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ABSTRACT 


The  design  of  a  system  devised  for  the  purpose  of  visually 
presenting  demographic  information  is  discussed  in  this  thesis. 

The  system  has  the  following  characteristics: 

1)  The  visual  presentation  is  in  the  form  of  maps,  each  of 
which  contains  a  map  outline  and  possibly  some  regional  boundaries, 
together  with  one  or  more  solids  in  line  drawing  form. 

2)  Each  solid  approximates  by  its  height  the  magnitude  of  a 
particular  statistic  at  that  geographic  location  that  the  solid's 
position  on  the  map  represents. 

3)  All  the  solids,  or  symbols,  on  a  map  are  height  variations 
on  one  basic  symbol  that  is  user -defined ,  but  convex,  with  no 
curves  or  dents. 

4)  A  three-dimensional  effect  is  conveyed  through  perspective 
projection  with  hidden  line  removal. 

5)  Map  sequences  are  produced  for  the  purpose  of  creating 
movies  that  show  the  fluctuation  of  statistics  with  time  at 
different  geographic  locations. 

6)  For  each  solid  the  user-provided  statistical  information 
must  be  in  terms  of  initial  and  terminal  solid  heights  for  a 
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given  length  of  movie  time.  For  any  particular  frame  the 
solid  height  is  obtained  by  linear  interpolation  between 
these  initial  and  terminal  values. 

7)  Maps  may  be  rotated  in  three-dimensional  space  as  the 
user  requires. 

The  aim  of  movie  generation  is  to  present  demographic 
material  in  the  form  of  a  time-lapse  study.  While  each  map 
provides  a  comparative  evaluation  at  different  locations, 
map  sequences  enable  comparative  evaluations  for  one  location 
at  different  times. 

The  bulk  of  this  thesis  deals  with  an  algorithm  for 
hidden  line  elimination.  In  order  that  the  movies  be  practical 
to  produce,  the  required  computer  time  must  be  minimized. 

Making  the  hidden  line  removal  process  generally  applicable 
and  yet  fast  was  the  demanding  objective  in  creating  the  system 
that  this  thesis  describes. 

The  system  description  is  implicitly  divided  into  two 
sections.  The  first  describes  how  user-provided  information 
is  used  to  create  the  original  maps  without  hidden  line 
removal.  And  the  second  details  the  methods  used  to  efficiently 
remove  hidden  lines.  Prior  to  this  description,  there  is  a  short 
discussion  of  general  mapping  techniques,  and  leading  off  after 
the  introduction,  a  motivational  analysis  of  this  pro ject^ that 
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is  a  whimsical  study  of  possible  application  areas. 
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CHAPTER  1 


The  Capabilities  of  the 
Dynamic  Data  Display  System 
as  Described  by  this  Thesis 
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This  thesis  describes  a  package  of  programs  that  was 
constructed  for  the  purpose  of  visually  presenting,  in  line 
drawing  form,  demographic  information.  Such  information 
provides  statistics  that  are  associated  with  graphical  areas 
or  locations.  The  system  enables  these  statistics  to  be 
visually  demonstrated  by  creating  maps,  with  particular 
locations  on  the  map  representing  the  geographic  locations 
to  which  the  demographic  information  to  be  processed  refers. 

The  purpose  of  such  a  map  is  to  somehow  reflect  statistical 
variations.  At  any  point  in  time,  the  magnitude  of  a  particular 
statistic  changes  from  one  geographic  location  to  another,  and 
a  map  can  be  composed  that  demonstrates  this  by  having,  at 
every  one  of  its  locations  that  corresponds  to  a  chosen 
location  within  a  geographic  area  under  study,  a  solid  or  symbol, 
that  reflects  the  magnitude  of  the  particular  statistic  by  its 
dimensions.  In  implementing  this  thesis,  only  symbol  height 
was  allowed  to  vary  in  accordance  with  statistical  fluctuation. 
Since  height  cannot  be  accurately  measured  visually,  it  is  a 
concept  of  relative  statistical  magnitude  that  is  conveyed. 
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In  order  that  visual  comparisons  between  symbols  be  as 

i 

accurate  as  possible,  every  symbol  is  constructed  by  altering 
the  height  of  one  particular  symbol  shape  that  is  user-defined. 
This  basic  symbol  must  be  convex,  with  no  dents  or  curved 
surfaces . 

To  clarify  in  the  reader's  mind  the  structure  of  a.  typical 
map,  consider  as  example  a  particular  study  that  has  Canadian 
population  as  the  statistic  under  consideration,  and  parallele¬ 
pipeds  as  the  required  symbols.  Imagine  a  traditional  map  of 
Canada  lying  on  a  table  top  with  matchboxes  standing  on  end  at 
those  locations  that  are  involved  in  the  population  study.  A 
sample  map  would  be  analogous  to  a  two-dimensional  representation 
of  this  with  all  but  the  outline  of  the  traditional  map  removed, 
and  with  those  matchbox  edges,  or  parts  of  edges,  that  cannot 
be'  seen  from  the  variable  position  in  3-space  from  which  the 
table  top  is  being  viewed,  also  removed.  Perspective  is 
maintained  in  the  two-dimensional  rendering. 

The  magnitude  of  a  statistic  not  only  varies  from  place 
to  place  within  a  geographic  area  at  any  one  time,  but  it 
also  changes  at  different  instances  in  time  for  any  one 
location.  The  programs  produce  maps  one  after  another  in 
such  a  way  that  changes  in  symbol  size  at  any  one  location 
on  the  map  reflect  statistical  fluctuation  at  the  corresponding 
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location  in  the  geographic  area. 
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In  following  a  pattern  of  going  from  the  general  to 
the  particular,  the  discussions  on  creating  the  maps  with 
hidden  lines  removed  are  delayed  by  a  discussion  of  general 
mapping  techniques,  with  some  indication  of  current  work 
involving  maps  and  the  computer,  and  also  delayed  by  a 
chapter,  to  which  the  reader  now  turns,  that  could  well  be 
called  a  flight  of  fancy.  And  to  the  reader  who  feels  that 
such  a  turn  of  fancy  has  no  place  in  an  endeavour  of  this 
nature,  I  ask  for  indulgence  on  the  grounds  that,  if  I  describe 
my  motives,  they  might  prove  contagious. 


CHAPTER  2 


Helping  Statistics  to  be  Heard  - 
In  Support  of  Computerized  Mapping 


"  'For  ever  shall  we  be  in  quest  of 
the  shores,  that  we  may  sing  and  be 
heard.  But  what  of  the  wave  that 
breaks  where  no  ear  shall  hear?  It 
is  the  unheard  in  us  that  nurses  our 
deeper  sorrow.  Yet  it  is  also  the 
unheard  which  carves  our  soul  to  form 
and  fashions  our  destiny. 1  " 


-  Kahlil  Gibran 
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In  March  of  1970,  a  conference  was  held  at  Evanston, 
Illinois,  to  introduce  to  the  world  of  computing  the  hopes 
of  geographers  and  cartographers  that  graphics  could  be 
harnessed  as  a  tool  for  their  ends.  The  feasibility  of 
such  a  tool  had  been  proven  by  a  Harvard -originated  package 
that  enabled  mapping  limited  to  gray  levels.  Significantly, 
the  conference  took  on  a  complexion  that  should  help  dis¬ 
courage  the  suspicions  that  many  seem  to  feel  toward  the 
scientist,  that  he  is  a  creature  dedicated  to  his  channelled 
purpose  and  oblivious  to  the  needs  and  sorrows  of  his  society. 
For  it  was  said  that  in  the  end  it  was  a  social  conscience 
that  was  heard  the  clearest,  bringing  to  light  the  potential 
value  of  computerized  mapping  to  those  social  problems  that 
could  be  visually  represented,  such  problems  as  pollution, 
traffic  noise,  and  unemployment.  Waves  of  statistics  presented 
by  word  of  mouth  tend  to  refract,  some  say,  by  the  densest  of 
all  mediums,  rather  than  produce  reflective  results,  but 
visually  presented  the  same  statistics  have  a  better  chance 
of  effecting  their  significance.  As  a  result,  computerized 
mapping  may  give  various  prognosticators  a  quiver  full  of 
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arrows  for  their  target  practice,  and  urban  planners  would 
derive  from  it  a  useful  tool  for  demonstrating  their  ideas. 

This  is  to  ignore  its  utility  to  many  other  organizations, 
such  as  police  departments  in  determining  deployment  of 
forces  to  various  parts  of  a  metropolis  as  crime  dictates 
the  need,  or  even  PTA  groups  for  using  childbirth  figures 
for  visualizing  that  blondes  really  do  have  more. 

It  should  not  be  thought  that  mapping  did  not  exist 
before  the  advent  of  the  computer.  What  it  contributes, 
however,  if  mapping  techniques  can  be  computerized,  is  speed 
in  analyzing  and  plotting  statistics.  It  becomes  practical 
to  plot  weekly,  even  daily,  statistical  changes,  assuming 
easy  collection  of  the  relevant  figures.  And  at  monthly 
intervals,  it  should  ultimately  be  possible  to  create  a  computer¬ 
generated  movie,  even  in  3-D,  that  would  help  provide  insight 
into  short-range  progressions.  What  a  computer  adds  is 
performance,  both  in  terms  of  speed  and  bulk  output.  It  would 
also  encourage  the  production  of  maps  depicting  data  on  less 
obvious  yet  socially  relevant  subjects,  such  as  the  number  of 
economically  active  females  over  age  15,  as  G.  M.  Gaites  of 
the  U.K. ,  Ministry  of  Housing  and  Local  Government,  has  done. 

It  might  be  said  that  data  is  today  being  collected  for  so 
many  different  subjects,  that  a  compiled  representation  of  all 
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our  statistics  would  effectively  reflect  what  could  be  called 

* 

the  quality  of  life.  This  quality,  as  shown  by  advanced 
mapping  techniques,  would  do  much  more  to  elicit  response 
than  facts  presented  by  digits  alone.  In  other  words,  mapping 
would  be  an  interface  between  the  cold  world  of  statistics  and 
the  general  social  conscience,  and  would  thus  considerably 
arouse  the  awareness  of  society  to  the  vicissitudes  of  the 
lives  that  make  it  up.  In  the  case  of  such  a  problem  as 
pollution,  though  it  is  carried  by  men,  it  is  contracted  by 
nature  and  thus  severely  effects  not  only  our  lives  but  those 
to  come.  Though  the  effects  of  pollution  have  been  slow  in 
materializing  to  the  proportions  we  know  today,  man's  reactions 
to  them  have  been  even  slower,  and  surely  the  responsibility 
of  each  person  is  not  only  to  realize  and  respond,  but  to  give 
others  the  time  to  react  also.  Here  computer  mapping  should 
play  an  important  role.  With  a  heightening  of  social  awareness, 
and  the  quality  of  life  thereby  in  closer  supervision,  responses 
to  its  ebb  and  flow  should  be  quicker,  and  we  may  not  call 
this  project  what  we  call  it  today,  a  mapping  for  survival. 

The  kind  of  mapping  that  would  probably  do  most  to  arouse 
social  consciousness  is  dynamic  mapping  -  showing  on  film  the 
changes  in  certain  statistics  through  a  time  period.  For 
example,  the  growth  in  population  over  a  span  of  years  in  some 
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part  of  the  world  could  be  well  demonstrated  in  this  manner. 

i 

Dynamic  mapping  would  be  conducive  to  the  extrapolation  of 
statistics  into  the  future,  and  this  representation  of  the 
problem  would  surely  speak  louder  than  words. 

The  production  of  films  to  show  changes  in  statistics 
could  be  very  valuable  as  an  educational  tool.  Visually 
presented  statistics  could  be  interpreted  to  show  changes  in 
social  attitudes,  economic  upheavals,  power  evolvements,  and 
general  social  evolution.  It  is  interesting  to  conjecture 
about  being  able  to  map  the  changes  in  the  responses  of 
society  to  polls,  and  thus  presenting,  on  screen,  trends  in 
thought,  in  feeling,  in  attitude.  The  study  and  understanding 
of  history  involves  as  much  the  study  of  its  socio-economic 
precursors  as  it  does  the  study  of  the  events  that  make  it. 

And  if  indeed  history  repeats  itself,  mapping  social  thought 
may  be  destined  to  serve  mankind  in  some  significant  way.  To 
realize  this  potential  role,  the  accessability  of  mapping 
means,  and  the  practicality  of  using  them,  both  in  terms  of  time 
and  money,  must  be  such  that  the  full  use  of  mapping  is 
encouraged.  Computers  can  rise  to  the  occasion. 

One  can  envision  a  mapping  package  that  could  satisfy 
just  about  any  application.  A  user  would  specify  the  outlines 
of  a  map,  the  method  of  statistical  representation,  whether  it 
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be  "band",  "field",  "spot"  (as  described  in  jis]  ),  3-D 

t 

parallelopiped ,  or  any  other;  he  would  also  specify  how  the 
statistics  should  be  scaled,  and  what  other  options  he  requires, 
such  as  entering  the  actual  numeric  form  of  the  statistics  on 
the  map,  entering  names,  or  in  the  case  of  a  movie,  perhaps 
entering  a  clock  or  calendar  to  represent  the  passing  of  time 
as  the  statistics  make  their  changes. 

Such  a  versatile  package,  though  the  ultimate  aim  of 
computerized  mapping,  is  not  the  subject  of  this  thesis.  My 
feelings  are  that  representing  the  larger  social  ills  is  of 
greatest  priority,  and  to  this  end  dynamic  mapping  in  three 
dimensions  is  the  most  effective.  This  form  of  representation 
is  not  the  most  appropriate  in  all  cases,  but  if  the  public 
is  going  to  stand  up  and  take  heed,  the  message  must  be  put 
forth  in  the  most  entertaining,  eye-catching,  and  if  science 
will  forgive  the  term,  artistic  way.  Imagine  the  world's 
food  supply  as  tall  parallelopipeds ,  full  and  stately  in  their 
imagined  sturdiness,  then  a  time  study  shows  them  shortening 
like  burning  fuses,  toppling  and  turning  in  deathly  throes, 
until,  extrapolating  into  the  next  century,  they  gradually 
diminish  to  the  stature  of  a  toy  Stonehenge  in  the  garden  of 
doll's  house,  dim  remnants  of  a  promising  past.  Or  population 
statistics  starting  out  as  uncluttered,  harmless  little 
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parallelopipeds ,  slowly  growing,  slowly  fattening,  strangling 
the  available  space,  threatening  each  other,  and  finally 
reaching  for  a  sky  that  cannot  possibly  hold  them,  ironically 
a  sky  leading  to  a  vast  and  empty  universe. 

We  can  only  conjecture  as  to  what  we  could  learn  if 
statistics  could  stand  up  and  present  themselves  so  that  their 
messages  could  be  heard.  Mapping  can  give  them  the  chance. 


That  is,  if  men  have  ears. 


CHAPTER  3 


The  Choice  of  Mapping  Technique 
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To  get  a  good  perspective  on  the  choice  of  the  mapping 
technique  that  is  the  subject  of  this  thesis,  a  brief 
description  of  the  cartographer's  art  is  in  order,  following 


the  lines  of  Fisher's  discussion 


For  our  purposes,  we  define  a  map  as  a  visually  readable 
display  of  quantitative  material  that  is  spatially  variable. 
The  method  of  representation  of  the  -material  is  the 
symbolism,  of  which  there  are  three  basic  types:  spot,  band, 
and  field. 

Spot  symbolism  consists  of  discreet,  compact  elements 
called  spots.  Usually  a  spot  is  centered  in  a  zone  that  it 
statistically  represents. 


Fig.  1  -  Spot  Symbolism 
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Band  symbolism  consists  of  elongated  and  somewhat  narrow 
bands  that  run  through  approximately  the  middle  of  the  regions 
to  which  they  apply. 


Fig.  2  -  Band  Symbolism 

With  field  symbolism,  the  symbols  are  relatively  widely 
spread  in  the  enclosing  region  so  that  they  form  a  field. 


Fig.  3  -  Field  Symbolism 
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Now  for  each  of  these  types  of  symbolism  there  are  three 

» 

techniques  for  quantitative  differentiation,  and  these  are 
extent,  tone,  and  countability. 

In  the  case  of  extent,  the  different  symbols  vary  in  some 
dimensional  aspect,  such  as  length,  width,  or  area.  For 
example,  with  spot  symbolism  the  greater  the  width,  area,  or 
volume  of  a  particular  spot,  the  greater  is  the  magnitude  of 
the  statistic  that  it  tries  to  represent.  For  bands  and  fields 
volume  is  not  appropriate,  although  height  and  proximity  may 
be  used. 

In  the  case  of  tone,  the  darkness  of  the  symbolism  varies 
directly  with  the  magnitude  of  the  statistic.  This  technique 
is  self-explanatory  in  all  three  cases  of  symbolism. 

With  the  method  of  countability,  the  number  of  elements 
within  the  symbolism  varies  directly  with  the  magnitude  of  the 
statistic.  It  is  necessary  that  the  number  and  arrangement  of 
the  elements  be  such  that  the  countability  is  visually  easily 
determined.  If  the  symbolism  is  via  spot,  then  a  particular 
spot  represents  a  value  that  is  directly  related  to  the  number 
of  elements  within  that  spot.  The  same  applies  to  bands  and 
fields . 

And  these  are,  essentially,  the  main  mapping  variations. 
For  an  extended  discussion  of  each,  refer  to  [is]  in  the 
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Each  variation  has  its  own  merits,  and  for  each  application 
certain  techniques  are  more  appropriate  than  others.  Ultimately 
a  mapping  package  may  be  created  that  would  enable  a  user  to 
chose  one  of  several  variations,  but  this  thesis  tackles  three- 
dimensional  spot  symbolism  in  association  with  quantitative 
representation  by  extent.  The  choice  was  partly  predicated  on 
the  challenge  of  the  hidden  line  problem,  which  would  only 
occur  with  the  use  of  three-dimensional  symbols.  It  was  also 
felt  that  the  greatest  visual  impact  would  be  obtained  by  such 
symbols,  although  it  should  be  clearly  stated  that  they  are 
useless  for  some  applications,  and  inappropriate  for  many  others. 
Yet  interest  in  the  general  field  of  computer  mapping  would  be 
best  aroused  by  as  flamboyant  an  initial  presentation  as  possible, 
without  making  it  unrealistically  hard  to  construct. 

And  this  it  is  not.  There  is  a  restriction  to  convex 
solids  that  ensures  this,  because  such  solids  have  no  curves, 
concavities,  or  irregularities  to  complicate  the  issue. 

Ultimately,  it  would  be  a  useful  achievement  to  be  able  to  work 
with  three-dimensional  depictive  shapes  as  symbols,  such  as 
humans,  animals,  cars,  etc.,  but  this  would  be  a  time-consuming 
pursuit  at  the  current  level  of  the  field  of  graphics. 

There  are  more  practical  reasons  for  chosing  a  spot-extent 
variation.  To  begin  with,  excluding  the  complication  of  solving 
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hidden  line  problems,  let  us  consider  the  advantages  and  dis¬ 
advantages  with  the  choice  of  extent  over  countability  and 
tone,  and  spot  over  band  and  field. 

First,  countability  requires  multiple  symbols  for  the 
zones  that  make  up  the  map,  while  extent  allows  for  one  symbol 
that  varies  in  size  with  the  value  of  the  statistic  that  it 
represents.  And  tone  is  also  more  complicated  to  implement  in 
that  not  only  the  outline  of  the  symbolism  is  under  consideration, 
but  also  all  the  interior  points.  If  the  mapping  format  is  made 
to  adhere  to  the  limitations  of  using  a  line  printer  as  the 
output  agent,  these  complications  may  not  be  so  important.  But 
for  paper  and  microfilm  plotters,  as  with  CRT's,  the  increase 
in  plotting  instructions  results  in  a  considerably  larger 
expense.  Since  the  creation  of  movies  was  the  ultimate  aim  of 
this  thesis,  extent  was  chosen  as  the  more  practical  method  of 
representing  statistical  variation. 

The  choice  of  spot  symbolism  was  also  based  partly  on  the 
fact  that  the  other  two  choices,  bands  and  fields,  tend  to  result 
in  a  greater  bulk  of  graphic  instructions.  This  is  because  the 
boundaries  of  bands  and  fields  are  determined  by  the  shape  of 
the  regions  they  represent,  and  map  regions  are  seldom,  if  ever, 
predictable  and  regular  in  outline.  In  the  case  of  spot  symbolism 
the  regions  are  irrelevant  except  to  locate  where  the  symbols 
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are  to  be  placed  -  the  shape  of  the  symbols  is  unchanged  from 

i 

region  to  region.  The  choice  of  symbol  determines  the  number 
of  graphic  instructions. 

Thus  the  spot-extent  variation  seems  to  have  the  best 
chance  of  success  when  considering  machine-oriented  advantages. 
However,  it  also  has  several  practical  disadvantages  that  are 
well-known  to  cartographers. 

The  first  concerns  the  variations  in  zone  sizes.  If 
the  zone  size  range  (which  is  the  difference  in  size  between 
the  largest  and  smallest  zones)  is  reasonably  large,  then  a 
particular  standard  symbol,  if  set  in  terms  of  the  smallest 
zone,  may  be  too  small  for  the  larger  zones.  And  if  the 
standard  is  set  for  the  largest  zone,  then  it  may  be  too  big 
for  the  smaller  zones.  Another  problem  that  is  linked  with 
zone  size  involves  visual  psychology.  Two  symbols  of  the 
same  dimensions  in  two  zones  of  greatly  differing  size  will  not 
look  identical  in  extent.  The  symbol  in  the  larger  zone  will 
seem  slightly  smaller  than  it  really  is,  while  the  reverse  is 
true  in  the  smaller  zone. 

Zone  size  is  also  important  when  there  is  a  large  variation 
in  the  values  of  the  statistics,  especially  in  a  situation  in 
which  a  time  study  is  involved.  Consider,  for  example,  using 
bars  to  plot  population  figures.  At  City  A  the  population 
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increases  extremely  rapidly  over  100  years,  while  at  City  B 
the  increase  is  moderate  and  the  final  figure  for  B  is  25 
times  smaller  than  for  A.  The  large  increase  at  A  has  to 
be  evident,  yet  the  final  bar  at  A  needs  to  be  within  size 
limits.  A  logarithmic  progression  will  not  give  good  visual 
accuracy;  yet  a  bar  at  B  25  times  smaller  than  a  bar  at  A  will 
mean  that  the  bar  at  B  will  look  insignificant.  If  area  is 
the  varying  factor  the  bar  at  A  must  be  reasonable  in  size 
relative  to  the  region  in  which  it  is  contained,  but  its  region 
may  be  much  smaller  than  the  region  containing  B.  Such  factors 
can  severely  complicate  the  mapping  problem,  but  though  they 
are  most  noticeable  in  the  spot-extent  mapping  variation,  they 
are  by  no  means  restricted  to  it. 

Another  issue  is  the  total  number  of  zones  in  the  map.  The 
greater  the  number  of  zones,  the  harder  it  becomes  to  visually 
distinguish  the  different  symbols.  In  genera],  the  three- 
dimensionality  of  symbols  will  lose  its  effectiveness  if 
there  are  too  many  zones. 

These  are,  then,  some  of  the  considerations  involved  in 
the  choice  of  mapping  technique.  In  the  final  analysis,  the 
choice  of  convex  solids,  changing  in  size  with  given  statistical 
values,  was  largely  based  on  the  intuition  that  such  a  system, 
within  the  context  of  film-making,  would  be  as  effective  as 
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possible  without  over-stepping  the  bounds  of  reasonable  simplicity. 

Others  have  combined  mapping  and  computers  in  different 
ways.  J.  D.  Jacobsen  {_12~|  has  created  a  system  that  manages 
to  retrieve  particular  information  associated  with  geographic 
areas,  given  a  data  base  that  represents,  in  terms  of  a  co¬ 
ordinate  system,  a  map  of  streets,  intersection  points,  zoning 
districts,  etc.,  about  which  there  is  a  geographic  fact 
information  that  may  include  population  statistics,  traffic 
density,  health  data*  and  other  demographic  information.  The 
user  has  the  ability  to  retrieve  statistical  data  from  a 
portion  of  a  map  within  the  two-dimensional  boundaries  of  a 
polygon  whose  construct  he  specifies.  Or  he  can  retrieve 
information  pertaining  to  some  particular  (X,  Y)  coordinate 
location.  The  important  characteristic  of  this  system  is  that 
it  is  designed  for  retrieving  geographically  oriented  data 
rather  than  a  system  for  representing  such  data.  As  such,  a 
very  obvious  application  area  is  in  the  design  of  an  urban 
management  information  system. 

SYMAP  is  a  program  designed  for  architects  and  planners 
that  was  developed  at  the  laboratory  for  Computer  Graphics  at 
Harvard  [13J  .  It  is  capable  of  producing  three  kinds  of 

maps:  (1)  chloropleth  or  conformant  maps  in  which  data  zones 

conform  to  particular  boundaries?  (2)  contour  maps,  based  on 
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contour  lines,  each  of  which  represents  a  uniform  value  through¬ 
out  its  length;  and  (3)  proximal  maps  based  on  proximity  to 
a  data  point.  A  cost  saving  feature  is  that  maps  are  output 
by  printer.  Characters  are  used  to  produce  tone,  with  the 
amount  of  over-printing  determining  the  degree  of  darkness. 

Tobler  (  [10  and  Qli|  )  has  been  working  on  models 
for  population  growth,  resulting  in  a  movie  that  simulates 
just  that  in  the  Detroit  region,  using  three-dimensional 
contour  maps.  At  the  University  of  Saskatchewan  work  is  being 
done  on  interactive  aids  for  compiling  maps  [l4]  .  To  my 

knowledge,  however,  very  little  work  has  been  done  in  creating 
maps  that  present  demographic  information  with  the  use  of 
convex  solids  with  hidden  lines  removed. 


CHAPTER  4 


Implementation  Tools 
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This  chapter  will  briefly  describe  the  tools  that  were 
used  to  implement  the  thesis.  The  essential  facility  was  an 
IBM  360/44,  operating  in  a  multi-programming  environment  with 
256  K,  and  having  attached  an  IBM  2250  Model  I  CRT  (with  key¬ 
board,  lightpen,  and  an  8  K  buffer) .  There  was  also  a  Calcomp 
835  microfilm  plotter  that  was  used  to  make  the  films. 

The  choice  of  software  proved  to  be  an  educational 
experience.  In  existence  at  the  time  that  the  thesis  was 
started  was  a  line  drawing  language  l/j  that  looked  quite 
suitable  for  the  dynamic  mapping  project.  This  language,  called 
AGILE,  was  intended  as  a  passive  aid  to  computerized  line  drawing, 
and  was  in  effect  a  set  of  primitive  graphic  functions  augmenting 
another  language  known  as  APL  [jL8~] 

There  are  several  central  characteristics  of  APL  that  makes 
it  so  suitable  as  a  base  for  AGILE'S  purposes.  To  begin 
with,  it  has  an  excellent  arithmetic  capability,  with  extensions 
of  the  usual  arithmetic  operations,  including  maximum  and  minimum 
value  capabilities,  random  number  generation,  etc.  Secondly, 
the  language  is  capable  of  working  with  matrices  as  operands. 

For  example,  A  +  B,  where  A  and  B  are  arrays  of  the  same  dimensions/ 
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results  in  an  array  also  of  the  same  dimensions  in  which 
each  element  is  the  sum  of  the  corresponding  elements  in  A 
and  B.  These  two  properties  among  others  were  thought  to  be, 
and  indeed  are,  useful  to  the  graphics  programmer  who  works 
with  pictures  formatted  in  the  classical  manner,  that  is,  with 
a  picture  represented  by  a  large  array  of  points  each  having 
associated  with  it  X,  Y  and  IC  values.  (The  IC  value  determines 
whether  a  visible  line  is  to  be  drawn  from  the  previous  point 
to  the  current  one) .  APL  has  good  capabilities  for  manipulating 
and  operating  on  arrays,  and  since  this  was  extended  by  AGILE 
to  operate  on  'picture'  arrays,  the  marriage  with  graphics 
seemed  likely  to  succeed. 

It  is  hard  to  determine  at  this  stage  whether  it  did.  The 

problems  that  arose  in  the  use  of  the  resultant  language  were 

caused  by  faulty  and  incomplete  implementation.  In  studying 

its  capabilities  it  looked  extremely  desirable,  and  made  light 

work  of  many  programming  requirements  where  other  methods  would 
# 

be  tedious  and  time  consuming.  But  in  operation,  AGILE  just 
gave  too  many  problems,  and  simply  did  not  fulfill  its  promises. 
So  it  was  dropped. 

Even  if  it  had  been  more  bug-free,  it  might  still  not  have 
been  suitable.  APL  has  the  characteristic  of  being  an  inter¬ 
preter,  having  been  initially  designed  as  a  conversational 
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language.  This  would  make  AGILE  extremely  slow  with  long 
compute-time  applications,  a  category  to  which  hidden-line 
removal  would  fall. 

At  any  rate,  FORTRAN  was  turned  to  as  the  basis  for 
dynamic  data  display.  The  probable  time  involved  in 
programming  in  ASSEMBLER  made  it  unattractive  (though  in  the 
long  run  it  would  be  most  efficient)  and  FORTRAN  was  chosen 
in  the  interests  of  time-saving  and  simplicity. 

A  program  that  simulates  graphic  output  on  a  line  printer 
was  used  extensively  during  the  testing  phases  to  minimize 


cost  and  turnaround  time. 


CHAPTER  5 


An  Overview 

of  the  Approach  Employed 


for  Data  Display 
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The  first  requirement  for  the  display  of  data  is  the  data 
itself,  with  some  way  of  organizing  and  inputting  it.  For 
this  thesis,  the  data  is  of  4  forms.  All  statistical  information 
will  be  represented  by  symbols  that  are  height  variations  of 
one  symbol  called  the  basic  symbol.  The  first  set  of  data 
describes  the  shape  of  the  convex  solid  acting  as  a  basic  symbol 
in  terms  of  the  coordinates  of  its  vertices  in  three-dimensional 
space  relative  to  the  origin.  Every  time  the  symbol  is  called 
upon  to  perform  a  representative  duty,  modifications  of  it 
will  include  a  change  of  height  and  of  orientation,  followed 
by  an  adjustment  of  the  coordinates  of  each  vertex  corresponding 
to  a  perspective  transformation  and  a  translation  that  will 
locate  the  new  symbol  at  its  appropriate  position.  It  is 
necessary  to  further  describe  the  basic  symbol  by  indicating 
the  number  of  faces  composing  the  symbol,  and  the  particular 
vertices  making  up  each  face,  and  yet  further  describe  which 
vertices  connected  by  straight  line  segments  construct  the 
symbol . 

Another  set  of  data  describes  the  points  making  up  the 
map  outline.  This  data  consists  of  the  coordinates  of  the 
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points  that  make  up  the  outline  together  with  an  indicator  for 
each  point  that  determines  whether  any  one  point  is  to  be 
joined  to  the  one  preceding  it. 

To  position  the  symbols  at  their  appropriate  locations 
requires  the  knowledge  of  the  coordinates  of  the  points  at  which 
the  symbols  are  to  be  placed.  Associated  with  each  such  point 
is  an  alphanumeric  identification,  so  in  effect  this  data 
corresponds  to  a  geographic  index  containing  a  list  of  place 
names  with  each  one  having  a  coordinate  code. 

The  final  set  of  data  is  made  up  of  the  statistics  that 
the  symbols  are  to  represent  through  their  size.  There  are 
two  questions  involved  in  the  aspect  of  change  of  symbol  size: 
what  function  is  to  be  used  in  relating  variation  in  statistics 
to  size  changes  (e.g.  linear,  logarithmic,  etc.),  and  how  a 
symbol  is  to  be  altered  to  visually  describe  such  changes 
(e.g.  modify  height,  volume,  etc) .  Both  questions  have  to 
be  resolved  in  the  light  of  the  following  considerations: 
minimize  symbol  cluttering,  help  simplify  the  viewer's  task 
of  visually  assessing  size  change,  and  not  to  omit  the  obvious, 
remember  that  the  frames,  all  of  which  are  the  same  fixed  size, 
have  to  be  capable  of  containing  the  largest  symbol  without 
making  the  smallest  unrecognizable. 
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For  the  purpose  of  implementation,  linear  height  variation 

« 

is  the  means  used  for  correlating  statistics  with  symbols,  and 
therefore  the  form  that  the  input  data  takes  is ,  rather  than 
raw  statistics,  starting  and  terminating  height  values  for  each 
symbol  for  a  given  length  of  movie.  Since  statistical  range 
varies  so  greatly  depending  on  the  study,  it  is  impractical  to 
allow  for  different  correlation  functions  between  statistics 
and  symbol  sizes  that  ensure  maximum  efficiency  for  all  studies. 

And  because  the  acceptable  convex  solids  are  numerous,  changing 
other  dimensional  aspects,  such  as  volume,  also  becomes  impractical. 
The  user  may,  however,  write  his  own  size  modification  program, 
particular  to  his  needs,  that  can  be  invoked  by  the  main  line  program. 
An  example  should  help  clarify  how  the  existing  display  system 
works.  It  is  based  on  the  following  configuration  of  co¬ 
ordinate  axes: 


-vv£  y 
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A  positive  rotation  about  an  axis  is  to  be  considered 

*  ■*  i 

as  a  counter-clockwise  rotation  when  looking  from  the  origin 
along  the  positive  part  of  that  axis. 

Consider  gathering  statistics  that  reflect  population 
growth  or  decline  over  the  last  60  years  in  two  mythical 
castles  that  shall  be  appropriately  called  Masterbegat  and 
Missbegat.  The  lord  and  his  lady  at  Masterbegat  have  been 
trying  for  30  years  to  have  female  offspring,  but,  alas,  their 
frantic  attempts  have  been  rewarded  by  no  less  than  30  males. 
At  Missbegat,  on  the  other  hand.  His  Lordship  has  matched 
stroke  for  stroke  his  neighbour's  potency  with  the  admirable 
result  of  a  YWCA  -  full  of  little  ladies,  also  30  in  number. 
Love  has  been  extremely  helpful  to  our  example  by  pairing  the 
existing  unmarried  eldest  of  each  brood  at  the  going  rate  of 
one  marriage  per  annum,  with  the  newlyweds  moving  to  shared 
accommodations  at  Togetherness  Castle.  This  matrimonial 
fait  accompli  began  one  year  after  the  30th  member  of  each 
brood  was  born,  and  during  the  whole  incredible  but  true 
events  of  the  entire  span  of  60  years  there  were  no  deaths, 
and  the  30th  member  of  each  family  was  definitely  the  last. 

If  each  castle  had  a  subservient  household  numbering  18,  then 
the  following  statistics  describe  the  population  changes: 
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YEAR 

MASTERBEGAT 

MISSBEGAT 

TOGETHERNESS 

0 

20 

20 

18 

1 

21 

21 

18 

2 

22 

22 

18 

• 

• 

n 

• 

20+n 

20+n 

18 

• 

30 

50 

50 

18 

31 

49 

49 

20 

32 

48 

48 

22 

i 

80-i 

80-i 

18+2x ( i-30) 

60 

20 

20 

78 

Now  suppose  that  a  60-second  movie  is  to  be  made  to 
visually  demonstrate  this  example. 

The  first  requirement  is  a  description  of  the  basic  symbol 
that  will  be  used.  In  keeping  with  the  use  of  parallelopipeds 
in  examples,  consider  a  cube  1  inch  x  1  inch  x  1  inch  as 
the  required  basic  symbol.  The  system  is  geared  to  accept  the 
coordinates  of  the  nodes  of  such  a  symbol  with  the  origin  at 
the  centre  of  the  symbol,  and  from  a  bird's  eye  vantage  point  - 
that  is,  looking  down  directly  at  the  top  of  the  symbol. 

Fig.  5.1  shows  how  the  cube  must  be  described. 
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The  Cube  as  a 
Basic  Symbol 


The  Point  of  View  and 
and  3-D  Orientation  Requir 
in  Describing  the  Basic  Symbol 


NODE  NUMBER  COORDINATES 

1  (-.5,  .5,  -.5) 

2  (-.5,  -.5,  -.5) 

3  (.5,  -.5,  -.5)  The  Description  of 

4  ( . 5 ,  .5,  -.5) 

5  (-.5,  . 5 ,  .5)  the  Basic  Symbol 

6  (-.5,  -.5,  .5) 

( • 5 ,  -.5,  .5) 

(.5,  .  5 ,  .5) 


Fig.  5.1 
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In  generating  map  symbols,  references  to  orientation 
and  height  are  made  with  respect  to  the  basic  symbol.  For 
example,  if  a  movie  is  to  start  by  showing  a  side  view  of  a 
particular  symbol,  it  is  generated  by  rotating  the  basic 
symbol  by  90  degrees  about  its  centre. 

The  map  outline  and  any  regional  boundaries  are  also 
initialized  using  a  top  view.  The  points  making  them  up 
should  initially  all  be  in  the  X-Y  plane.  The  locations  at 
which  the  symbols  are  to  be  placed  should  also  be  initially 
provided  set  into  the  X-Y  plane.  Using  the  IBM  835  Microfilm 
plotter  frame  dimensions  of  20  inches  x  15  inches,  let  the 

coordinates  for  the  castles  in  our  example  be  (in  inches) : 


(5,  5,  0) 
(10,  10,  0) 
(16,  6,  0) 


Masterbegat 

Missbegat 

Togetherness 


Now  suppose  that  the  movie  is  to  be  made  showing  the 
60  years  of  the  population  study  in  60  seconds.  Furthermore, 
the  initial  frame  is  to  show  the  parallelopipeds  in  silhouette 
or  side  view,  that  is,  with  each  rotated  by  +90  degrees  about 
the  X-axis  from  the  original  position  established  by  the 
basic  symbol.  Each  symbol  is  to  be  rotated  during  the  first 
30  seconds  in  two  directions  -  from  the  initial  position, 
about  the  X-axis,  to  a  position  +45  degrees  from  the  basic  symbol 
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orientation  with  respect  to  the  same  axis,  and  also  from  the 
initial  position,  about  the  Y-axis,  to  a  position  -45  degrees 
from  the  basic  symbol  orientation  with  respect  to  the  Y-axis. 

In  the  subsequent  30  seconds  the  symbols  are  to  return  to  their 
starting  orientations.  Fig.  5.2  shows  the  angle  changes. 


Time  =  0 

Side  View 
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Time  =  60  secs. 

Side  View 


Fig.  5.2 


At  this  point  it  is  appropriate  to  discuss  how  rotations 
are  handled.  All  elements  on  a  map  -  each  symbol  and  the  map 
outline  -  are  rotated  separately  by  the  same  amount.  A  symbol 
that  is  to  represent  some  demographic  information  is  created 
by  adjusting  in  linear  fashion  the  height  of  the  basic  symbol, 
rotating  it  about  its  centre,  and  placing  it  at  the  appropriate 
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location  on  the  map.  But  the  coordinates  of  this  location  must 
also  be  adjusted  according  to  the  rotation.  As  the  plane  of 
the  map  outline  shifts,  points  that  are  on  this  plane  must 
shift  with  it.  The  map  outline  is  rotated  about  the  centre  of 
the  frame,  and  the  locations  have  to  be  moved  in  a  corresponding 
manner  in  3-D  space.  However,  if  this  were  the  only  adjustment 
for  coordinates  of  symbol  locations,  then  a  side  view  (as 
required  in  the  example)  would  show  the  map  outline  as  a  solid 
line  half  way  up  the  frame  and  parallel  to  the  X-axis.  All 
symbols  would  be  standing  on  this  line  and  have  only  half  the 
frame  to  contain  them.  Adjustments  are  therefore  made  to  lower 
the  map  outline  and  symbols.  The  amount  of  adjustment  is  such 
that  for  a  full  side  view  the  outline  is  a  line  parallel  to 
and  one  inch  above  the  X-axis.  Fig.  5.3  demonstrates  the 
adjustment . 


Side  View  With 
Adjustment 


Side  View  Without 
Adjustment 


Fig.  5.3 
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Now  we  can  return  to  the  particular  example  at  hand. 

Each  parallelopiped  has  to  be  given  a  starting  and  ending 
height  for  each  of  the  two  parts  of  the  movie.  The  maximum 
height  that  any  parallelopiped  may  achieve  should  still  allow 
that  parallelopiped  to  fit  within  the  frame  dimensions.  It  is 
easiest  to  consider  a  full  side  view.  If  the  height  of  a  frame 
is  15  inches,  and  if  in  side  view  the  base  of  a  symbol  is  one 
inch  above  the  X-axis,  then  no  parallelopiped  may  attain  a 
height  greater  than  14  inches.  Let  us  arbitrarily  limit  the 
height  to  12  inches.  This  height  then  corresponds  to  the 
height  of  the  largest  possible  parallelopiped,  which  in  our 
example  will  be  at  Togetherness  castle  in  the  60th  year  when 
the  population  is  at  a  maximum  of  78.  A  scale  is  then  calculated 
at  12/78  =  0.154  inches  per  person  -  that  is,  each  individual 

adds  this  height  to  a  parallelopiped.  For  the  first  30  seconds, 
then,  the  height  changes  are  (in  inches)  : 


CASTLE  STARTING  HEIGHT  ENDING  HEIGHT 


Masterbegat 

.154 

X 

20  =  3.08 

.154 

X 

50  =  7.70 

Missbegat 

.154 

X 

20  =  3.08 

.154 

X 

50  =  7.70 

Togetherness 

.154 

X 

18  »  2.77 

.154 

X 

18  =  2.77 
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And  for  the  last  30  seconds: 


CASTLE 

STARTING  HEIGHT 

ENDING  HEIGHT 

Masterbegat 

.154 

x  49  =  7.55 

.154  x  20  =  3.08 

Missbegat 

.154 

x  49  =  7.55 

.154  x  20  =  3.08 

Togetherness 

.154 

x  20  =  3.08 

.154  x  78  =  12.01 

We  can  distinguish  two  parts  to  the  data  needed  for  creating 
a  movie:  data  that  applies  to  the  movie  as  a  whole,  and  data 
that  applies  to  sub-sections  of  the  movie.  In  the  former 
category  there  is  the  description  of  the  basic  symbol,  the  des¬ 
cription  of  the  map  outline  and  the  regional  boundaries,  the 
number  of  symbols  to  be  generated  per  frame,  the  data  required 
to  locate  the  symbols  on  the  map,  and  the  coordinates  of  the 
viewpoint. 

In  the  second  category,  there  are  two  batches  of  input 
required  for  our  example,  one  for  each  half  of  the  movie.  In 
summary,  the  information  that  is  required  for  the  generation 
of  the  first  30  seconds  is: 

1)  The  length  of  movie  time:  30  seconds. 

2)  The  starting  and  ending  heights  of  each  symbol: 


SYMBOL  # 

STARTING 

HEIGHT 

ENDING 

HEIGHT 

1 

3.08 

inches 

7  .70 

inches 

2 

3.08 

inches 

7.70 

inches 

3 

2.77 

inches 

2.77 

inches 
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The  starting  map  angle  (as  already  described) : 

i 

a)  +90  degrees  about  the  X-axis  from  the  initialized 
orientation. 

b)  0  degrees  about  the  Y-axis. 

c)  0  degrees  about  the  Z-axis. 

The  ending  map  angle  (as  already  described) : 

a)  +45  degrees  about  the  X-axis  from  the  initialized 
orientation. 

b)  -45  degrees  about  the  Y-axis. 

c)  0  degrees  about  the  Z-axis. 


The  following  steps  will  be  taken  by  the  system  to  create 
this  part  of  the  movie: 


1)  Find  the  number  of  frames  to  be  created. 

No.  of  frames  =  30  x  24  =  720 

2)  '  Find  the  height  change  per  frame  for  each  symbol. 

e.g.  at  Masterbegat  this  rate  is: 

(7.70  -  3.08)/720  =  .0064  inches/frame. 

3)  Find  the  angle  changes  per  frame. 

(45  -  90)/720  =  -.0625  degrees/frame  about  the  X-axis, 
(-45  -  0)/720  =  -.0625  degr ees/frame  about  the  Y-axis, 

0  degrees/frame  about  the  Z-axis. 

4)  Create  a  frame  at  a  time  till  720  frames  are  made.  For  each 
frame,  the  following  steps  are  necessary,  using  the  300th  frame 
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in  examples: 

a)  Find  the  map  angles  for  the  current  frame.  With 
reference  to  the  initialized  orientation,  for  the  300th  frame, 
these  are: 

90  +  (-.0625)  x  300  =  71.25  degrees  about  the  X-axis, 

0  +  (-.0625)  x  300  -  18.75  degrees  about  the  Y-axis, 

0  degrees  about  the  Z-axis 

b)  Rotate  the  outline  by  the  appropriate  angles  and 
adjust  the  height  to  lower  the  outline  if  there  is  rotation 
about  the  X-axis.  Then  apply  perspective  projection.  The 
method  for  doing  this  is  discussed  later. 

c)  As  in  b) ,  adjust  the  coordinates  of  the  locations 
at  which  the  symbols  are  to  be  placed. 

Consider  the  symbol  for  Togetherness  Castle.  Its  original 
coordinates  are  (16.00,  6.00,  0.00).  Rotation  about  the  centre 
of  the  frame  at  (10.00,  7.50,  0.00)  by  the  angles  found  in  a) 
results  in  the  coordinates  (16.15,  7.02,  0.57).  Since  rotations 
in  3-space  are  not  commutative,  it  should  be  noted  that  the 
order  of  rotation  is  first  about  the  X-axis,  then  the  Y-axis, 
and  finally  the  Z-axis. 

The  coordinates  of  each  location  must  also  be  lowered 
to  provide  symbols  with  maximum  space  in  which  to  grow,  as 
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has  already  been  discussed.  The  Y  coordinate  is  reduced  by 
6.5  x  (1.0  -  cosine  of  the  angle  of  rotation  about  the  X-axis). 
For  a  side  view,  the  map  outline  will  be  a  line  one  inch  above 
and  parallel  to  the  X-axis,  since  the  frame  midpoint  has  co¬ 
ordinates  (10.0,  7.5).  For  Togetherness,  the  location  co¬ 
ordinates  now  becomes  (16.15,  1.92,  0.57) . 

d)  Each  symbol  is  generated  in  turn.  To  demonstrate 
the  process,  consider  the  symbol  for  Masterbegat.  The  steps 
in  creating  it  are: 

i)  Create  a  new  symbol  by  modifing  the  height  of 
the  basic  symbol.  The  height  change  is  made  symmetric  about 
the  X-Y  plane.  For  the  300th  frame  the  height  is  3.08  x  300  x 
4.62/720  =  5.0  inches.  The  coordinates  of  the  nodes  then  are: 


NUMBERING 

SCHEME 

NODE 

NUMBER 

COORDINATES 

1 

'v 

/ 

9 

1 

(-.5,  .5,  -2.5) 

’i  s 

2 

(-.5,  -.5,  -2.5) 

i 

1 

3 

(.5,  -.5,  -2.5) 

i 

1 

l 

4 

(.5,  .5,  -2.5) 

I 

it 

1 

5 

(-.5,  .5,  2.5) 

2 

-  -  / 

6 

(-.5,  -.5,  2.5) 

2 

2 

7 

(.5,  -.5,  2.5) 

8 

(.5,  .5,  2.5) 

Note  that 

the 

magnitudes  of 

the  Z  coordinates 

of  the  basic 

symbol  nodes 

are  irrelevant. 

Only  their  signs 

are  necessary. 

ii) 

The  symbol  is 

now  rotated  about 

the  origin  (still 
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at  its  centre  )  by  the  appropriate  angles. 

iii)  The  important  task  now  arises  of  placing  the 
symbol  at  the  appropriate  position  on  the  map. 

Assuming  no  perspective  projection,  this  is  done  in  two 
steps.  The  first  is  to  add  to  the  Y  coordinate  of  the  existing 
symbol  location,  the  absolute  value  of  one-half  the  Y  component 
of  the  symbol  height.  This  is  to  move  the  symbol  upwards  by 
an  amount  dependent  on  the  angle  of  rotation  about  the  X-axis 
until,  with  a  full  side  view,  the  symbol  stands  with  its  base 
on  the  outline  (see  Fig.  5.3) .  An  adjustment  of  this  kind  is 
necessary  because  otherwise,  for  one  fixed  orientation  (excluding 
a  top  view)  a  symbol  would  appear  to  grow  symmetrically  about 
its  centre  when  it  should  only  be  rising  at  its  top.  The  second 
step  is  to  add  to  the  coordinates  of  each  node  of  the  symbol, 
the  corresponding  location  coordinates  in  order  to  place  the 
symbol  at  the  required  map  location. 

With  perspective  projection,  the  second  step  is  slightly 
altered  so  that  it  modifies  the  X  and  Y  coordinates  of  each  node, 
according  to  certain  projection  formulae,  while  leaving  the  Z 
coordinate  uneffected.  A  line  of  sight  is  considered  from  the 
viewpoint  to  the  "middle  of  the  X  -  Y  coordinate  plane.  In  order 
to  implement  the  projection,  therefore,  the  symbols  must  be 
provided  with  their  node  coordinates  relative  to  this  middle  point. 
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This  is  done  by  modifying  the  location  coordinates  to  have  this  point 
as  its  origin,  and  only  then  adding  to  the  node  coordinates 
the  corresponding  location  coordinates.  A  picture  plane  is 
now  placed  in  3-space  perpendicular  to  the  line  of  sight  and 
in  such  a  location  that  all  nodes  of  all  symbols  will  always 
be  on  the  opposite  side  of  the  plane  as  the  viewpoint.  The 
symbol  is  projected  onto  this  plane  and  is  then  translated 
to  restore  the  origin  to  its  prior  location. 

This  method  is  also  applied  to  the  map  outline  and 
regional  boundaries  (in  step  4b) .  The  projection  formulae 
are  derived  in  Appendix  C. 

When  all  the  symbols  for  a  frame  are  generated,  hidden 
lines  are  removed  from  them,  excluding  the  outline  and  regional 
boundaries.  This  process  is  described  later. 

When  the  first  half  of  the  movie  has  been  created,  the  data 
for  the  second  half  must  be  read  and  processed  in  the  same 
manner . 

The  example  provided  is  about  as  simple  as  possible,  while 

■ 

also  being  practical.  The  use  of  more  sophisticated  methods 
(e.g.  logarithmic  size  variation,  statistical  magnitude  related 
to  volume,  etc.)  will  require  the  user  to  program  size  change 
as  a  subroutine.  Since  each  symbol  has  its  own  peculiarities, 

it  is  not  practical  to  construct  a  programming  package  to  allow 

' 

for  all  contingencies.  However,  greater  choice  of  size  variation 
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methods  has  to  be  considered  a  worthwhile  extension  to  this 

i 

thesis.  Flexibility  in  this  department  was  in  effect  sacrificed 
for  a  good  working  solution  to  the  hidden  line  removal  problem 
and  flexibility  in  user-defined  symbolism. 


CHAPTER  6 


How  Pictures  are  Structured 
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This  chapter  will  describe  the  data  structure  used  to 
define  and  manipulate  symbols,  and  then  the  format  employed  to 
represent  a  composite  of  symbols,  that  together  with  a  map 
outline,  make  up  a  picture. 

The  basic  function  of  the  dynamic  mapping  set  of  programs 
is  to  accurately  project  onto  a  two-dimensional  surface  a 
three-dimensional  world.  When  we  view  a  particular  three- 
dimensional  environment,  there  is  almost  always  some  element 
in  it  that  is  in  part,  if  not  in  whole,  obscured.  It  is  the 
most  barren  of  environments  in  which  this  is  not  the  case.  In 
projecting  this  environment  on  a  two-dimensional  surface,  in 
order  to  maintain  realism,  methods  have  to  be  developed  to 
accurately  reproduce  those  parts  that  are  visually  perceptible 
in  3-space  and  eliminate  the  others,  as  seen  from  a  particular 
viewpoint.  This  has  been  aptly  named  a  problem,  specifically 
the  "hidden  line  problem"  or  "hidden  surface  problem" . 

In  order  to  describe  the  data  structure,  it  is  necessary 
to  have  adequate  terminology  to  describe  components  of  the  3-D 


and  2-D  worlds. 
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Definitions  and  Nomenclature 


1.  A  O-dimensional  polytope  1  1"0  is  a  point. 

2.  An  n-dimensional  polytope  1  1^  in  Euclidean  space  is  a 

finite  set  of  polytopes  with  the  following  properties, 

-tl  1  : 


a)  Every  Tlrv-x  which  is  an  element  of  a  I  1  -rv-i  i  T 
belongs  to  exactly  one  other 


*Y\. 


-i  £ 


TV.  • 


b)  No  subset  of  "J  forms  a  polytope 


Tv. 


of  its 


own. 

3.  A  line  segment  is  a  1-dimensional  polytope  or  a  pair  of 
points . 

4.  A  polygon  is  a  2-dimensional  polytope  or  a  set  of  line 
segments  such  that  the  beginning  point  of  one  is  the  end  point 
of  exactly  one  other. 

5.  A  vertex  is  a  O-dimensional  polytope  that  makes  up  either 
end  of  any  line  segment  of  a  polygon. 

6.  A  polyhedron  is  a  3-dimensional  polytope  made  up  of  a  set 
of  polygons  joined  at  common  line  segments. 

7.  An  edge  is  that  part  of  a  polyhedron  that  projects  onto 
the  viewing  plane  as  a  line  segment. 

8.  A  node  is  an  endpoint  of  an  edge  of  a  polyhedron.  (The 
next  chapter  will  redefine  node  into  2  classes  -  original  and 


unoriginal) 
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9.  A  face  is  that  part  of  a  polyhedron  that  projects  onto 
the  viewing  plane  as  a  polygon. 

i 

Thus,  on  the  viewing  plane,  nodes  are  projected  into 
vertices,  edges  into  line  segments,  and  faces  into  polygons. 

We  will  now  delve  into  the  structural  representation  of 
polyhedra  by  the  user.  Only  convex  polyhedra  are  considered. 
Definitions 

The  interior  INT  (P)  of  a  polytope  P  is  that  finite  portion 
of  Euclidean  space  -that  the  polytope  encloses,  while  the 
boundary  Bd  (P)  of  a  polytope  P  is  the  union  of  all  the  interiors 

of  the  polytopes'  various  polyhedra  or  k-cells^  0  ^  ^  ^  Y\.  —  \  : 

~\ — r  i  fc-L  -t— r'u. 

Bd  (  I  Ire  )  =  .U  U  INT  (  |  |  •  ) 

1-0  L 

A  convex  polytope  P  is  a  polytope  TT-w  ,  Tv  >  0  such  that 
for  any  points  XI,  X2  £  INT  (P)  U  Bd  (P) ,  the  interior 
of  the  line  segment  X1-X2  lies  entirely  in  the  interior  of  P. 

Only  convex  polytopes  are  allowable  as  symbols.  Intuitively, 
then,  we  only  accept  symbols  with  no  curves  and  no  dents. 

It  is  the  user  who  must  define  the  basic  symbol  that  is 
to  be  used  to  represent  all  statistical  information.  There  are 
3  aspects  to  the  user's  description.  The  first  provides  a 
description  of  the  shape  of  the  symbol  by  means  of  the  coordinates 
of  the  nodes  of  the  symbol  with  reference  to  the  origin,  as 
detailed  in  the  previous  chapter.  The  second  aspect  is  information 
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that  indicates  the  structure  of  the  faces  of  the  symbol,  namely 
the  number  of  faces  comprising  each  symbol  and  the  nodes  that 
in  turn  make  up  each  face.  Finally,  detail  is  provided  to 
enable  the  system  to  draw  between  pairs  of  vertices  the  lines 
that  make  up  the  symbol  outline.  The  second  and  third  aspects 
of  symbol  definition  are  actually  much  the  same,  but  the 
programming  is  made  less  complex  by  the  inclusion  of  both. 

Examining  each  aspect  of  symbol  definition  in  turn,  it 
is  evident  that  the  nodes  of  a  symbol  must  be  ordered.  Though 
explicit  node  numbering  is  not  needed,  implicitly  the  first 
three  coordinates  are  considered  to  belong  to  node  1,  the  next 
three  to  node  2,  and  so  on.  The  method  of  describing  a  basic 
symbol,  and  how  it  is  manipulated,  is  discussed  in  the  previous 
chapter  and  will  not  be  repeated  here.  However,  for  illustrative 
purposes,  the  diagram  of  the  cube  as  a  basic  symbol  is  exhumed. 


Fig.  6.2  -  A  Cube  as  a  Basic  Symbol 
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The  second  description  required  for  the  basic  symbol 
involves  the  faces  of  the  symbol.  Referring  to  Pig.  6.2,  it 
is  evident  that  there  are  6  faces,  with  the  nodes  making  them 
up  following  this  pattern: 


FACE  NUMBER 
(Arbitrary) 

1 

2 

3 

4 

5 

6 


NODES  OF 
FACE 

1,  2,  3,  4 

2,  6,  7,  3 
5,  6,  7,  8 

1,  5,  8,  4 

2,  6,  5,  1 

3,  7,  8,  4 


Note  that  .the  vertices  of  the  polygon  formed  by  projecting 
a  face  on  the  viewing  plane  have  to  be  ordered  clockwise  or 
counter-clockwise  about  a  ray  from  the  viewpoint  to  an  interior 
point  of  the  polygon.  The  ordering  may  be  different  for  different 
faces,  however. 

With  some  symbols,  not  all  the  faces  contain  the  same 

number  of  nodes.  The  maximum  number  of  nodes  per  face  is  required, 

> 

and  the  description  of  all  other  faces  must  allow  for  this  maximum. 
A  0  is  used  for  padding  where  necessary,  since  a  Oth  node  does 
not  exist.  For  example,  consider  the  symbol  of  Fig.  6.3. 
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I 


Fig.  6.3 

For  Fig.  6.3  there  are  5  faces  whose  node  configuration 
description  would  be  placed  in  a  5  x  4  matrix  as  follows: 

12  5  0 

2  3  4  5 

13  4  0 

15  4  0 

12  3  0 

To  complete  the  symbol  description,  again  referring  to 
Fig.  6.2,  the  system  must  know  which  nodes  to  join  by  straight 
lines.  Where  a  0  appears  in  this  connectivity  description,  the 
light  beam  is  turned  off  when  moved  to  the  succeeding  vertex. 

If  two  zeroes  appear,  the  description  is  assumed  ended.  Either 
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one  of  the  two  following  strings  of  node  numbers  could  serve 
the  desired  purpose  for  the  cube  in  Fig.  6.2: 

123702658067840341500 

or 

123415678502603704800 

For  the  first  case,  the  steps  that  would  be  followed  in  creating 
a  parallelopiped  would  be  as  in  Fig.  6.4. 


Group  1 

vertices  1,  2,  3,  7 


\ 


t 

1  7 


Group  3 

vertices  6,  7,  8,  4 


Group  2 
vertices  2,  6 


5,  8 


Group  4 


vertices  3,  4,  1,  5 


/  f 


71 


Z./ 


t'  ~Y3 


5: 

/ 


I  ! 


7* 


/ 


/ 


/ 


1/ 


 .V 


1/ 


7 


Fig.  6.4 
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The  user  should  be  aware  that  the  most  efficient  format 
for  the  connectivity  description  is  the  one  that  minimizes  the 
path  length  in  tracing  out  the  symbol,  including  the  relocations 
required  after  a  0  in  the  description. 

When  the  coordinates  of  the  nodes  of  all  the  symbols  that 
are  to  make  up  the  final  picture  have  been  calculated  along 
the  lines  described  in  the  previous  chapter,  they  are  placed  in 
the  columns  of  a  3  by  N  array  INAR.  If  J  is  the  number  of  nodes 
in  the  basic  symbol  to  be  employed  in  a  particular  movie  sequence 
and  if  there  are  K  symbols,  each  a  copy  of  the  basic  symbol,  then 
N  =  J  x  K,  with  the  nodes  of  the  first  symbol  in  the  first  J 
columns,  the  nodes  of  the  second  in  the  succeeding  J  columns, 
and  so  on.  The  reason  that  coordinates  are  placed  in  columns, 
not  rows,  is  that  the  coordinates  of  any  one  point  will  be  in 
contiguous  core  locations,  and  this  simplifies  subroutine 
invocation  when  the  coordinates  of  a  point  are  to  be  passed  as 
parameters . 

The  internal  representation  of  a  picture  to  be  displayed 
also  follows  this  pattern,  with  the  exception  that  the  third 
coordinate  of  a  point  is  replaced  by  a  value,  called  an  IC  value, 
that  indicates  whether  the  point  is  to  be  joined  by  a  straight 
line  to  the  one  preceding  it  in  the  picture  array.  Such  an 


array  looks  like: 
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XI  X2  X3 . XN 

i 

Y1  Y2  Y3  . .  .YN 

IC1  IC2  IC3 . . . . . ICN 


This  picture  format  is  one  that  is  used  by  a  graphics 
package  at  the  University  of  Toronto  devised  by  L.  Mezei 
called  SPARTA  jjL6]  .  It  was  important  that  the  dynamic 
mapping  programs  should  be  compatible  with  SPARTA  in  the 
event  that  some  of  the  picture  manipulation  properties  of 
SPARTA  were  to  be  exploited.  A  slight  variation  on  the 
array  described  makes  it  completely  acceptable  to  SPARTA 
routines.  This  variation  is  the  addition  of  a  column  at  the 
start  of  the  array  with  3  elements  that  have  to  be  initialized 
with  certain  data.  The  first  element  must  contain  the  maximum 
length  of  the  array  (N,  where  the  array  dimensions  are  3  x  N) . 
The  second  element  contains  the  value  1.0  if  the  initialization 
procedure  has  been  executed,  and  the  third  contains  a  number 
that  indicates  the  existing  extent  of  the  picture  that  the 
array  represents,  including  the  control  information  being 
described.  This  extent  refers  to  the  number  of  active  columns 
and  must  not  be  greater  than  N. 

Now  that  we  have  a  picture,  hidden  lines  have  to  be  removed 


from  it. 


CHAPTER  7 


Hidden  Line  Removal 
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Solutions  to  the  problem  of  hidden  line  removal  have 
evolved  very  little  from  those  propounded  when  computer 
graphics  was  first  introduced  to  the  mechanized  world.  I 
am  basing  this  view  on  two  facts.  First,  current  methods 
cannot  be  applied  satisfactorily  to  any  arbitrary  object 
in  3-D  space.  And  second  is  the  pragmatic  consideration 
that  removing  hidden  lines  from  all  but  the  simplest  objects 
is  a  time  consuming  and  tedious  endeavour,  given  current 
algorithms.  The  algorithm  used  for  the  dynamic  mapping 
project  is  one  that  restricts  the  objects  to  convex  solids, 
but  strongly  attempts  to  minimize  hidden  line  removal  time. 

There  are  two  very  distinct  groups  of  methods  for 
removing  hidden  lines.  The  first  is  more  accurately  called 
hidden  surface  elimination. 

This  first  group  of  methods  aims  at  creating  pictures 
in  chiaroscuro,  in  which  3-D  effects  are  generated  through 
the  use  of  shading.  Devices  are  not  currently  capable  of 
generating  subtle  changes  in  tone,  so  that  the  best  that  can 
be  attempted  is  the  simulation  of  the  half-tone  process  of 
printing.  Several  important  and  worthwhile  attempts  have 
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resulted,  and  the  reader  is  referred  to  Q3]  ,  JVJ  ,  jjjJ 


The  Wylie  et  al  version  |J3J  approximates  all  surfaces 
with  small  but  finite  triangles.  An  object  so  divided  and 
its  perspective  projection  on  the  viewing  plane  are  scanned 
by  a  ray  extending  from  a  given  viewpoint.  The  scan  ray 
passes  over  the  entire  picture  by  scanning  horizontally, 
dropping  by  a  vertical  increment  AY,  scanning  horizontally 
again,  and  so  on.  Every  time  the  ray  drops  by  AY  a  Y -occupied 
table  is  updated  to  reflect  those  triangles  that  the  scan  ray  will 
potentially  intersect  in  its  new  horizontal  trajectory.  An 
X-occupied  table  keeps  track  of  those  triangles  that  are 
actually  being  intersected  by  the  ray  at  the  current  (X,Y) 
location.  This  table  is  updated  at  every  /\X  if  any  triangles 
are  exited  or  entered  in  making  the  horizontal  increment.  If 
more  than  one  triangle  is  intersected  by  a  ray,  a  hidden  parts 
calculation  determines  which  face  having  that  portion  of  its 
surface  bounded  by  the  edges  corresponding  to  the  line  segments 
making  up  the  triangle  sides,  is  closest  to  the  viewpoint  along 
the  ray  from  the  viewpoint  to  the  point  under  test.  A  visible 
table  records  the  results  for  a  particular  scan  line.  It  is 
then  necessary  to  calculate  the  intensity  per  point  of  the 
scan  line.  The  intensities  at  the  vertices  of  the  component 
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triangles  are  calculated  before  the  entire  scanning  process, 

i 

using  fundamental  laws  of  light  reflection,  assuming  the 
viewpoint  as  the  source.  The  intensity  of  a  point  within  a 
triangle  is  obtained  through  interpolation  of  the  intensities 
at  the  three  vertices.  Thus  the  intensity  associated  with  a 
point  along  a  scan  line  is  calculated  using  the  vertex 
intensities  of  the  triangle  that  the  visible  table  claims 
contains  that  point.  The  list  of  intensities  for  the 
scan  line  is  sent  to  a  peripheral  device  for  eventual  display. 

This  is  a  representative  method  of  creating  half-tone 
pictures,  so  no  other  will  be  described.  The  point  to  notice 
about  the  method  is  that  it  examines  the  entire  surface  of 
an  object  rather  than  the  edges  that  make  it  up,  which  is 
the  case  with  the  second  group  of  hidden  line  removal  methods. 
This  second  group  produces  line  drawings  rather  than  chiaroscuro 
pictures . 

While  half-tone  pictures  require  the  examination  of  every 
point  on  every  face  in  a  picture,  line  drawings  are  produced 
by  reducing  edges  into  visible  and  invisible  portions,  and 
this  requires  finding  points  of  intersection  of  edges  with 
planes,  and  line  segments  with  line  segments  (where  a  line 
segment  is  the  projection  of  an  edge  onto  the  viewing  plane) . 

Basically,  the  method  of  removing  hidden  lines  in  a  line 
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drawing  is  an  exercise  in  geometry.  In  its  most 
nrimitive  form,  every  point  in  3-space  that  belongs  to  any 
edge  is  checked  for  obscurity  by  all  faces,  using  3-D 
geometry.  This  brute  force  method  is  just  too  time  consuming, 
so  most  algorithms  use  shortcuts  that  help  reduce  the  amount 
of  hidden  point  checking.  These  shortcuts  attempt  to  locate 
areas  of  line  segment  intersection  and  treat  them  separately 
from  line  segment  and  edge  endpoints.  In  general,  no  curved 
faces  are  allowed.  Some  algorithms  simplify  matters  by  only 


removing  hidden  lines  from  isolated  objects 


7 


Most 


often  the  objects  themselves  are  restricted  to  disallow  any 
indentations  (convex  objects) .  A  few  algorithms  employ 
unorthodox  methods,  eg.  |j3  J  ,  but  the  bulk  seem  to  use  the 
brute  force  approach  with  time  saving  modifications  that  can 
be  applied  to  a  certain  restricted  class  of  pictures.  As 
I  have  said,  the  algorithm  to  be  described  here  is  applicable 
to  convex  polyhedra  that  may  intersect  in  3 -space  with  their 


own  kind. 
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THE  HIDDEN  LINE  ALGORITHM 

The  following  definitions  and  nomenclature  are  required 
before  describing  the  algorithm. 

Definitions  and  Nomenclature 

A  picture  of  objects  that  lie  in  Euclidean  space  (3-space) 
is  that  polygon  on  a  given  viewing  plane  that  only  contains  all 
the  polygons  formed  by  projecting  the  faces  of  the  objects  on 
the  same  viewing  plane. 

An  object  will  be  said  to  be  in  a  picture  if  every  polygon 
formed  by  projecting  the  faces  of  the  object  onto  the  viewing 
plane  is  within  the  polygon  that  defines  the  picture  on  the 
same  viewing  plane. 

An  original  node  of  a  picture  is  any  point  in  3-space  that 
is  an  endpoint  of  any  edge  of  any  object  in  the  picture.  The 
set  of  all  original  nodes  of  a  picture  will  be  designated  by 


N,  where  N  is  the  set 


Note 


that  all  subscripts  reference  a  set  member.  Nj  £  N,  but 
NI  is  any  node  that  may  or  may  not  be  a  member  of  N.  All 
nomenclature  referring  to  vertices,  line  segments,  edges,  polygons, 
and  faces,  follows  the  same  pattern. 
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An  unoriginal  node  of  a  picture  is  any  node  that  is 

i 

not  original  but  that  lies  on  any  edge  of  any  object  in  the 
picture. 

The  picture  space,  PS ,  of  a  picture  is  the  set  of  all 
original  and  unoriginal  nodes  of  the  picture. 

An  edge  joining  any  two  original  nodes  Nj  ,  Nj,  will 
be  designated  Nj  -  Nj.  If  the  nodes  may  be  unoriginal,  then 
they  are  designated  NI  and  NJ  and  the  edge  becomes  NI  -  NJ. 

For  any  Npr  £  N  ,  VK  is  the  vertex  formed  by  projecting 
NK  on  the  viewing  plane.  For  NK  £  PS,  the  vertex  is  VK. 

The  line  segment  formed  by  projecting  an  edge  Nj  -  Nj  on 
the  viewing  plane,  Nj  £  N,  Nj  £  N,  will  be  designated 
Vj.  -  Vj.  If  the  nodes  are  unoriginal,  the  line  segment  is 
VI -VJ . 

IT  is  the  set  of  all  faces  of  the  objects  in  a  picture. 
The  members  of  F  are  designated  by  F]_  ,  F2  ,.... . ...F^. 

The  polygon  formed  by  projecting  FK  on  the  viewing 
plane  will  be  POLY^. 

S_  will  refer  to  the  set  of  all  solids  or  symbols  that 
are  involved  in  making  up  a  picture.  Subscripts  again  denote 
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members  of  this  set. 

F(SR)  refers  to  the  set  of  all  faces  of  SR,  SR  £.  S. 

A  face  Fj  belongs  to  SK  iff  Fj  £  F(Sk)  ,  SK  £  S. 

N(Fk)  is  the  set  of  all  original  and  unoriginal  nodes  of 

fK'  fk  e  F- 

N(Sk)  is  the  set  of  all  original  and  unoriginal  nodes  of 
SK'  SK  £•  s- 

A  node  NI  belongs  to  SR  iff  NI  £,  N(SR)  ,  NI  i  PS,  SK  ^  S. 

fj(Sk)  is  the  set  of  nodes  making  up  the  Jth  face  of 
symbol  SR,  Fj  £  F,  SR  £  S. 

An  edge  pierces  a  face  FR,  FR  £_  F,  if  the  intersection  of 
the  edge  and  face  in  3-space,  when  projected  on  the  viewing 
plane,  lies  inside  the  boundaries  of  POLYR. 

PLANE R  is  that  plane  in  3-space  that  continues  face  FR, 

fk  £  F- 

VP  is  the  point  in  3-space  that  serves  as  the  viewpoint. 

A  point  in  3-space  will  be  said  to  be  hidden  from  VP  by 
a  face  FR,  FR  £,  F,  if  the  point  lies  on  the  opposite  side 
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of  Fk  as  VP,  and  if  it  lies  within  the  boundaries  of  POLYK. 

i 

A  visible  edge  is  one  seen  from  VP  in  its  entirety. 

If  an  edge  NI  -  NJ,  NI  £.  PS,  NJ  £,  PS ,  does  not  pierce 
a  face  FK,  fk£F,  but  if  VI  -  VJ  intersects  POLYp^  at  VL,  then 
the  node  formed  by  projecting  VL  on  NT  -  NJ  will  be  designated 
by  P (VL,  NI  -  NJ) . 

A  face  Fk,  Fk  £.  F  ,  is  on  view  if  no  part  of  its  surface 
is  hidden  from  VP  by  any  face  Fj  where  Fp^  £  F(Sk)  /  Fj  £.  F(Sk)  ' 
SK  £  s« 

The  term  list  will  be  used  to  describe  a  structure  that 
looks  like  a  matrix,  but  which  can  have  multiple  information 
as  an  element.  For  example,  an  element  may  contain  a  node, 
although  a  node  requires  three  coordinates  to  describe  it.  We 
will  arbitarily  say  that  a  row  of  a  list  is  an  entry,  and  an 
element  of  an  entry  will  refer  to  an  element  in  the  list  that 
belongs  to  that  entry.  An  entry  may  be  added  to  or  taken  away 
from  a  list  at  will. 
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At  this  point,  we  review  the  information  that  is  given 
about  a  picture  for  which  hidden  lines  are  to  be  removed. 


(1)  The  coordinates  of  all  the  original 
Consider  a  picture  that  contains  5 


the  original  nodes  are: 


Ni  N2 
v _ 


e  • 


>n8n9, 


■Nl6 


nodes  of  the  picture, 
parallelopipeds .  Then 


s5 


N  =  [n],  ,  N2  ,  N3  . N  40*1 

s  =  [sx  ,  s2  ,  s3  ,  s4  ,  s 5 "y 

Associated  given  information  therefore  includes  the  number  of 
original  nodes  per  symbol. 

(2)  A  vector  describing  how  the  original  nodes  of  any  one 
symbol  are  to  be  joined  to  form  its  edges. 

If  for  the  example  in  (1)  this  vector  is: 
123702658067840341500 
then  the  edges  of  the  picture  are: 
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EDGES  OF  Sx  EDGES  OF  S2  EDGES  OF  S5 


Nl 

-  n2 

N  9  -  N10 

n33 

-  n34 

n2 

-  n3 

N10  “  N11 

n34 

-  n35 

N3 

-  N7 

Mu  -  N15  . 

n35 

"  N39 

N 

2 

N6 

N10  -  N14 

N34 

7  N38 

N,  „  - 
12  9 


N  9  -  N14 


S-^  looks  like: 


W, 

r 


IS/* 


U' 


H  5  L 


/ 


VrN. 


^6 


n 


n7- 


Fig.  7.2 
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(3)  A  matrix  that  determines  the  inter-relationship  between 
faces  and  nodes  for  any  one  symbol. 

If  the  matrix  looks  like: 


then, 

are: 


1 

2 

3  4 

2 

6 

7  3 

5 

6 

7  8 

1 

5 

8  4 

1 

2 

6  5 

3 

7 

8  4 

again 

us  ing 

the 

example  of 

(1) 

Fi(Si) 

— 

NX, 

n2. 

n3,  n4 

P2(Sl) 

1 

= 

n2. 

N6, 

N7,  N3 

• 

ft 

f8<si) 

= 

n3. 

• 

• 

Ny , 

n8,  n4 

f1<S2) 

= 

n9. 

N10 

,  Nlx,  1 

w12 

f2(s2) 

% 

1 

N10' 

N14;  N15' 

Mu 

1 

f8(s5) 

=3 

n35# 

9 

N39,  n40. 

n36 

the  faces  of  the  picture 


Note  that  the  original  nodes  of  any  face  FK,  FK<^F, 
must  be  ordered  in  some  circular  fashion,  with  a  node  recorded 
only  once. 

The  hidden  line  algorithm  will  now  be  described  in  several 
phases . 
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PHASE  I 

i 

This  phase  handles  4  initialization  procedures. 

A.  Find  'Common  Boxes  1 

If  two  symbols  Sj  and  Sj,  Sj  £,  S,  Sj  £.  S,  intersect  in 
3-space,  then  their  common  box  is  that  'box'  in  3-space  whose 
faces  are  all  parallel  to  the  coordinate  planes  and  which 
has  within  or  on  its  face  boundaries  that  part  of  Sj  that 
intersects  in  3-space  with  S j.  Fig.  7.3  and  Fig.  7.4  demonstrate 
this  concept.  There  may  also  be  a  common  box  when  the  symbols 
do  not  intersect  in  3-space.  This  occurs  if  there  is  an  inter¬ 
section  of  the  projections  of  Sj  and  Sj  on  the  viewing  plane. 

In  this  case,  the  polygon  formed  by  projecting  onto  the  viewing 
plane  the  front  face  of  the  common  box  will  contain  the  area 
in  common  between  the  projections  of  Sj  and  Sj.  Fig.  7.5 
provides  an  example  of  this  kind  of  common  box. 
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Fig.  7.3 


A 

v\€w(:>oi*j'r 


Fig.  7.4 


V  \  (:  w  ^  i 
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Y  2 


V  \  6  u)  P  O  \ 


Fig 


7.5 
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Consider  a  common  box  CB.  The  information  required  to 

i 

specify  the  dimensions  of  the  common  box  is:  the  smallest 
and  largest  X  coordinate  of  its  original  nodes,  CBMINX  and 
CBMAXX  respectively,  the  smallest  and  largest  Y  coordinates, 
CBMINY  and  CBMAXY  respectively,  and  the  smallest  and  largest 
Z  coordinates,  CBMINZ  and  CBMAXZ  respectively. 


y 


(C0>nt*x  CO-mAv^) 


C  C^mInx,  cftntMV,  C0,m  Ax-tj 


i  (^CCbrlNy  f  t.  0  r  AxY,  ./ 


(cGnfixx,  C6hAxy;  tGpflv?) 


Ccf^HAxv^  c(ihA*y, 


(cGrftxx^  tdhrNy,  CCt^Avz) 


(.cG^Axx  t  to,  mIm  ^ 


y  cd,t-AX^v;  (.ftmiA'-'a) 


lt 


Fig.  7.6  -  Dimensions  of  a  Common  Box 
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Definitions 

A  node  NI  will  be  said  to  be  hidden  by  a  common  box  if 
it  is  hidden  from  VP  by  the  front  face  of  the  box,  or  if  VI 
lies  on  a  line  segment  that  is  a  boundary  of  the  polygon 
formed  by  projecting  the  front  face  of  the  common  box  on  the 
viewing  plane,  and  also  lies  behind  the  plane  continuing 
the  front  face. 

An  edge  will  be  said  to  intersect  a  common  box  if  any 
point  on  the  edge  is  hidden  by  the  common  box. 

The  procedures  for  determining  whether  a  point  is  hidden 
by  a  common  box,  and  whether  an  edge  crosses  a  common  box, 
do  not  merit  inclusion.  They  are  trivial  enough. 

Finding  the  common  box  of  two  symbols  Sj  and  Sj,  Sj  £  S, 
Sj  £.  S,  is  also  a  simple  matter.  That  'box'  with  faces 
parallel  to  the  coordinate  planes  is  found  which  has  within 
or  on  its  boundaries  every  N(Si) .  Such  a  'box'  is  called  the 
enclosing  box  of  Sj.  Say  that  the  maximum  X  coordinate  of 
its  nodes  is  SIMAXX,  the  maximum  Y  SIMAXY,  the  maximum  Z 
SIMAXZ,  the  minimum  X  SIMINX,  the  minimum  Y  SIMINY,  and 
the  minimum  Z  SIMINZ.  If  the  corresponding  dimensional 
requirements  for  Sj  are:  SJMAXX,  SJMAXY,  SJMAXZ,  SJMINX, 
SJMINY,  SJMINZ,  then  we  can  say  the  following  about  the 
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common  box  of  S-j-  and  Sjt 

a)  it  does  not  exist  if  any  of  the  following  hold  true 

i)  SIMAXX  <  SJMINX 

ii)  SIMINX  >  SJMAXX 

iii)  SIMAXY  <  SJKINY 

iv)  SIMINY  >  SJMAXY 

(Refer  to  Fig.  7.7) 


1/ 


Fig.  7.7 
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b)  If  no  condition  in  a)  holds,  then  the  common  box 

i 

specifications  are: 


CBMAXX 

= 

smaller  of 

SIMAXX 

and 

SJMAXX 

CBMAXY 

= 

smaller  of 

SIMAXY 

and 

SJMAXY 

CBMAXZ 

= 

smaller  of 

S IMAXZ 

and 

SJMAXZ 

CBMINX 

= 

larger  of 

SIMINX 

and 

SJMINX 

CBMINY 

= 

larger  of 

SIMINY 

and 

SJMINY 

CBMINZ 

= 

larger  of 

SIMINZ 

and 

SJMINZ 

(An  example 

is 

shown  in  Fig.  7.8) 

Y 


Sr 


4 


fixV  < 


A vy  <  SjhAx  v 


sa  m  Ax  V 


V 


Fig.  7.8 
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It  is  to  be  noted  that  a  common  box  does  not  always 
map  particularly  accurately  the  areas  of  interaction  between 
two  symbols.  For  example,  refer  to  Fig.  7.9  in  which  the  solids 
and  the  viewer  are  being  looked  at  from  an  aerial  vantage  point. 


Fig.  7.9 

Looking  down  at  2  parallelopipeds  and  a  viewer. 


The  degree  of  interaction  between  Sj  and  Sj  does  not  approach 
being  commensurate  with  the  common  box  dimensions.  Edge  NK  -  Nl 
of  Sj  is  considered  to  intersect  the  common  box  even  though  it 
has  nothing  in  common  with  Sj.  This  will  cause  some  needless 
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calculations.  But  the  advantages  to  having  the  common  box 
with  its  faces  parallel  to  the  coordinate  axes  outweighs 
extra  calculations  of  this  nature. 

For  any  two  symbols  Sj  and  Sj,  Sj  £  S,  Sj  £  S,  there 
is  an  entry  in  a  common  box  table  containing  the  box  speci¬ 
fications.  Zeroes  indicate  no  common  box.  The  table  saves 
considerable  time  in  two  aspects  of  hidden  line  determination: 
1)  If  a  node  NT  £  N(Sk)  ,  £  S,  is  to  be  checked  for 

obscurity  from  VP  by  an  on  view  face  of  SL,  SL  £  S,  then  we 
can  assume  no  possibility  of  NI  being  hidden  if  (a)  there  is 
no  common  box  between  SK  and  SL  ,  or  (b)  NI  is  not  hidden  by 
the  common  box,  if  it  exists.  A  point  that  is  on  the  boundary 
of  a  common  box  cannot  be  ignored  for  obscurity.  Referring 
to  Fig.  7.10,  node  N  is  hidden  even  though  it  lies  on  a 
side  face  of  the  box. 

) 


Fig.  7.10 
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2)  If  a  line  segment  VI  -  VJ  projected  from  an  edge  NI  -  NJ, 

i 

where  NI  £  N(Sk)  ,  NJ  £  N(SK)  ,  is  to  be  checked  for  an  inter¬ 
section  with  a  line  segment  VF  -  VG  projected  from  an  edge 
NF  -  NG ,  where  NF  £  N(SL) ,  NG  £  N(SL) ,  such  that  SR  £  S ,  SL  £  S, 
then  there  cannot  be  any  intersection  if  any  of  the  following 
hold  true: 

(a)  there  is  no  common  box  between  SK  and  SL. 

(b)  NI  -  NJ  does  not  intersect  the  common  box,  if  it  exists. 

(c)  NF  -  NG  does  not  intersect  the  common  box,  if  it  exists. 

B .  Find  On-View  Faces 

Each  allowable  symbol  will  have  one  or  more  faces  that 
are  on  view,  and  one  or  more  faces  that  are  entirely  hidden. 

Since  all  the  faces  of  a  picture  are  implicitly  numbered,  a 

vector  is  established  which  contains  in  the  Ith  element  a  1 

if  the  Ith  face  is  on  view  and  a  0  otherwise. 

To  discover  whether  a  face  of  a  symbol  SK,  SK  £  S,  is 
on  view,  find  a  point  in  3-space  that  is  an  interior  point 
of  S„.  This  point  is  given  an  X  coordinate  that  is  an 
average  of  the  X  coordinates  of  the  first  and  third  members  of 
N(Sk) ,  a  Y  coordinate  that  is  an  average  of  the  first  and  third 
members  of  N(SK)  and  a  Z  coordinate  that  is  an  average  of  the 
first  and  third  members  of  N(SR) .  A  face  of  SK  is  on  view 
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if  it  hides  the  interior  point  from  VP. 


C .  For  Every  On-View  Face,  Calculate  The  Values  That 

Determine  The  Equation  Of  The  Plane  in  Which  It  Lies 


The  equation  of  a  plane  is  given  by  AX  +  BY  +  CZ  =  D. 
Given  3  nodes  on  a  plane  whose  coordinates  are  (XI,  Yl,  Zl) , 


(X2 ,  Y2 , 

Z2)  , 

and 

(X3,  Y3,  Z3) ,  the  values  of  A, 

B, 

c. 

are: 

* 

Yl 

Zl 

1 

Zl 

XI 

1 

A  = 

Y2 

Z2 

1 

B  = 

Z2 

X2 

1 

Y3 

Z3 

1 

, 

Z3 

X3 

1 

XI 

Yl 

1 

XI 

Yl 

Zl 

C  = 

X2 

Y2 

1 

D  = 

X2 

Y2 

Z2 

X3 

Y3 

1 

X3 

Y3 

Z3 

For  every  on-view  face,  these  values  are  obtained  and  placed 
in  a  matrix. 

D .  Find  On-View  Edges 

This  is  a  simple  process  that  lists,  for  every  symbol, 
those  edges  of  a  symbol  not  hidden  by  any  part  of  itself. 
These  edges  are  just  those  that  belong  to  the  on-view  faces. 
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E .  Flag  Nodes  Belonging  To  Off-View  Faces 

A  vector  H  is  established  whose  length  corresponds  to 

the  number  of  original  nodes  in  the  picture  being  processed.  The  Ith 

element  of  this  -vector  will  eventually  contain  0  if  the  Ith 

node  is  hidden,  or  a  positive  integer  that  acts  as  a  code 

to  index  the  face  that  hides  the  node.  However,  if  the  node 

is  hidden  because  it  belongs  to  an  off-view  face,  then,  at 

this  stage,  a  -1  is  placed  in  the  appropriate  position  to 

indicate  that  the  node  may  be  ignored  in  future  calculations. 
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PHASE  II 

i 

In  this  phase  each  original  node  N j  £  N(Sj<;) /  £.  S, 

Nj  not  already  found  hidden  in  E  of  Phase  I,  is  checked 

for  obscurity  by  any  on-view  face  Fj  such  that  Fj  ^  F(Sk) • 

A  general  procedure  is  now  described,  that  can  be  used  to 
perform  this  function. 

The  General  Method  For  Finding  Whether  A  Node  NI  Is  Hidden 

By  A  Face  FK 

The  node  is  hidden  if  two  conditions  hold: 

a)  the  node  is  behind  PLANE^  as  seen  from  VP. 

b)  the  vertex  VI  lies  within  the  boundaries  of  POLY^. 

To  find  out  whether  a)  is  true,  get  the  value  obtained 
by  substituting  the  coordinates  of  NI  in  the  equation  of 
PLANE pr  (as  determined  in  Phase  I).  If  this  value  is  0,  the 
node  is  on  PLANE^  and  condition  a)  does  not  hold.  Otherwise, 
get  another  value  by  substituting  the  coordinates  of  VP  in 
the  plane  equation.  If  the  two  values  are  of  opposite  sign, 
the  node  is  behind  PLANEK, otherwise  condition  a)  does  not  hold 
and  the  node  is  visible. 

To  find  out  whether  condition  b)  holds,  it  is  first 
necessary  to  find  the  coordinates  of  any  point  that  lies  within 
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the  boundaries  of  POLYK.  The  X  and  Y  coordinates  obtained 
by  averaging  the  X  and  Y  coordinates  respectively  of  the  first 
third  vertices  of  POLYK  belong  to  an  internal  point.  If  the 
line  joining  VI  to  the  internal  point  intersects  any  boundary 
of  POLY^,  then  VI  is  outside  POLYK,  otherwise  it  lies  inside. 
This  method  was  chosen  because  the  process,  as  defined  in  the 
appendix,  that  determines  whether  two  line  segments  intersect, 
is  quick  and  painless. 

Returning  to  the  function  of  Phase  II,  if  N-j-  is  found 
to  be  hidden  by  an  on-view  face  FK'  the  code  that  refers  to 
Fk  is  placed  in  the  appropriate  element  of  the  vector  H 
introduced  in  E  of  Phase  I.  Fig.  7.11  shows  the  entries  in 
H  for  a  particular  two  symbol  picture. 

Phase  II  makes  the  first  use  of  the  common  box  concept. 
For  a  face  Fj,  Fj  £.  F(Sp)  ,  Fj  is  only  checked  for  obscuring 
the  node  Nj  under  consideration,  Nj  £_  N(Sk)  ,  if  there  is  a 
common  box  between  Sk  and  Sl,  and  if  N^.  is  hidden  by  the 
common  box  when  it  exists. 


and 
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NODE  NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


VALUES  IN  H 

0 

0 

0 

0 

-1 

7 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

0 
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PHASE  III 

In  this  phase  every  on-view  edge  belonging  to  every  symbol 
in  the  picture  is  worked  on  separately.  We  will  consider  a 
representative  edge  Nj  -  Nj,  whose  corresponding  members  in  H 
are  Hj ,  Hj,  and  where  Nj  -  Nj  belongs  to  SK '  SK  £  S* 

Three  possibilities  arise  for  on-view  edges: 

a)  Nj  and  Nj  are  both  hidden,  i.e.  Hj  >  0,  Hj  7  0. 

b)  One  node  is  hidden,  the  other  visible, 

i.e.  Hj  >  0  and  Hj  =  0,  or 

H-j.  =  0  and  Hj 

c)  Both  nodes  are  visible,  i.e.  Hj  .  =  Hj  =  0. 

Each  will  be  discussed  in  turn. 

A.  Both  Nodes  Hidden 

•  If  Ht  =  H T  both  nodes  are  hidden  by  the  same  face 
-L  U 

and  the  edge  is  totally  hidden. 

We  now  consider  either  node,  Nj  for  example.  Let  the 
polygon  formed  by  projecting  Hj  on  the  viewing  plane  be 
POLYHj ,  and  let  the  plane  continuing  POLYH-j-  be  PLANEHj .  The 
edge  is  now  processed  as  follows.  First  find  whether  Nj  -  Nj 
intersects  Hj  in  3-space.  We  do  this  in  two  steps, 
a)  Find  if  Nj  -  Nj  intersects  PLANEHj.  Substitute  N^  into 
the  equation  of  PLANEHj  getting  VALl,  and  substitute  Nj  into 
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the  same  equation,  getting  VAL2  .  Then: 

1)  If  VALl  and  VAL2  have  opposite  signs,  Nj  -  Nj 
intersects  PLANEHj. 

2)  If  VALl  and  VAL2  have  like  signs,  there  is  no 
intersection. 

3)  Since  a  node  is  visible  when  it  lies  on  a  plane, 
and  since  both  Nj  and  Nj  are  hidden,  then  neither  VALl  nor 
VAL2  will  be  0. 

b)  If  Nj  -  Nj  intersects  PLANEH-j-  then  3-D  geometry  is  used 
to  find  the  coordinates  of  the  point  of  intersection  (see 
Appendix  B) .  It  is  then  necessary  to  discover  whether  the 
projection  of  this  intersection  on  the  viewing  plane  is  within 
POLYHj  ,  using  the  method  demonstrated  in  Phase  II.  If  it  is 
within  POLYH-j-  ,  then  Nj  -  Nj  intersects  the  face  and  this 
intersection  point  replaces  N-j-  as  one  endpoint  to  form  a  new 
edge  which  will  be  called  NI  -  Nj.  (See  Fig.  7.13)  . 

If  an  intersection  in  3-space  is  not  found7  V-j-  -  Vj  will 
intersect  a.  line  segment  forming  a  boundary  of  the  projected 
polygon  of  .  If  this  intersection  is  VL,  a  new  endpoint, 

NI,  that  correspondents  to  P(VL,  Nj  -  Nj) ,  replaces  Nj. 

(see  Fig.  7 . 14) . 


y 
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N  T 


A' 


Fig.  7.13 


VI  f  PCs-j.  m  T 


Fig.  7.14 
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The  new  segment  NI  -  Nj  is  now  considered.  It  is  now 
necessary  to  find  whether  NI  is  hidden  from  VP,  and  this  is 
simply  a  matter  of  repeating  Phase  II  with  NI  only,  using 
the  same  available  shortcuts  that  the  common  box  concept 
allows.  If  NI  is  hidden,  there  is  equivalently  a  new  edge 
whose  endpoints  are  hidden,  and  a  return  is  made  to  A  of 
this  phase,  substituting  NI  for  N^and  the  new  hiding  plane 
for  Hj .  In  effect,  this  procedure  takes  that  part  of  the 
edge  starting  at  Nj  and  traces  it  till  it  emerges  into  sight 
at  some  NI .  At  this  point  there  is  an  edge  NI  -  Nj  with  one 
endpoint  hidden,  and  transfer  is  made  to  that  logic  which 
handles  such  edges. 

B .  One  Node  Hidden 

The  procedure  for  processing  an  edge  with  one  endpoint 
hidden  and  the  other  visible  is  identical  to  that  described 
in  A.  It  is  simply  a  question  of  tracing  the  edge  from  NT 
toward  NI  till  it  emerges  into  sight  at  some  NJ.  Now  there 
are  two  visible  endpoints  and  a  branch  is  made  to  the  logic 
that  handles  just  such  edges. 

C .  Both  Nodes  Visible 

Consider  an  on-view  edge  NI  -  NJ,  NI  and  NJ  belonging 
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to  Spr#  SK  £  S.  It  is  necessary  to  process  this  supposedly 
totally  visible  edge  by  first  finding  all  the  points  of 
intersection  of  VI  -  VJ  with  other  line  segments.  The  process 
is  as  follows.  Every  member  SL  of  S  is  taken  in  turn,  SL  ^  S^. 
If  SL  and  SK  do  not  have  a  common  box,  the  current  is 
ignored.  If  there  is  a  common  box,  SL  is  again  ignored  if 
NI  -  NJ  does  not  intersect  it.  Otherwise,  each  on-view  edge' 
of  SL  is  taken  in  turn.  If  any  such  edge  does  not  intersect 
the  common  box,  it  is  ignored.  Otherwise,  an  intersection  is 
found  between  VT  -  VJ  and  the  line  segment  formed  by  projecting 
the  edge  of  SL  onto  the  viewing  plane,  if  such  an  intersection 
exists.  Any  intersection  on  the  viewing  plane  is  placed  as 
an  entry  in  a  list  we  shall  call  NINT,  and  the  on-view  face 
or  faces  to  which  the  edge  belongs  is  recorded  as  an  entry 
in  another  list,  called  for  argument's  sake,  FACELIST.  Since 
an  edge  could  potentially  belong  to  two  on-view  faces,  there 
are  two  elements  per  entry.  If  there  is  only  one  face  to  be 
recorded,  one  element  is  made  to  contain  a  0  to  indicate  the 
situation.  Thus  the  Ith  entry  of  FACELIST  refers  to  the  Ith 
entry  of  NINT. 

When  all  intersections  have  been  recorded,  they  are 
sorted  by  increasing  distance  from  VI  and  replaced  in  NINT. 

Any  change  in  entry  order  for  NINT  is  applied  identically 


86 


to  FACELIST. 

i 

Definitions 

We  will  say  that  a  line  segment  interacts  with  a  polygon 
if  any  point  on  the  line  segment  lies  inside  or  on  the 
boundaries  of  the  polygon. 

A  segment  VI  -  VJ  enters  a  polygon  at  VK  if  VK  is  the 
first  point  at  which  VI  -  VJ  interacts  with  the  polygon  in 
tracing  the  path  from  VI  to  VJ.  VK  will  be  called  the  point 
of  entry. 

A  segment  VI  -  VJ  exits  a  polygon  at  VL  if  VL  is  the 
last  point  at  which  VI  -  VJ  interacts  with  the  polygon  in 
tracing  the  path  from  VI  to  VJ.  VL  will  be  called  the  point 
of  exit. 


VL 

• 

Pgl.V  * 

-  VT 


VI  —  VJ  enters  POLY^  at  VK.  VK  is  the  point  of  entry. 
VI  -  VJ  exits  POLYk  at  VL.  VL  is  the  point  of  exit. 


Fig.  7.15 
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A  modified  version  of  FACELIST  is  now  placed  in  another 
list  FI NT.  The  Ith  entry  of  FINT  is  made  to  reflect  the  faces 
whose  polygons  have  been  entered  but  not  exited  by  the  line 
segment  joining  NINTj-  to  NINTI+^  (the  length  of  the  list  FINT 
is  one  entry  less  than  the  length  of  NINT) .  Since,  in  effect, 
the  first  occurrence  of  a  polygon  in  the  Ith  entry  of  FACELIST 
heralds  a  point  of  entry  into  that  polygon  at  NINTj,  and  the 
second  occurrence  of  a  polygon  in  the  Jth  entry  of  FACELIST 
heralds  a  point  of  exit  from  that  polygon  at  NINTj,  the 
transition  from  FACELIST  to  FINT  is  extremely  easy.  Fig.  7.16 
provides  an  example. 


V  r '  -  88  - 

C  icrvae  (3la^£ 


NINT 

FACELIST 

FINT 

Comments 

PT1 

POLYl 

POLYl 

POLYl 

is  entered  at  PT1 

PT2 

POLY2 

POLYl , 

POLY  2 

POLYl 

POLY  2 

has  not  been  exite 
is  entered  at  PT2 

PT3 

POLYl,  POLY 3 

POLY2 , 

POLY3 

POLYl  is  exited  at  PT3 
and  is  therefore  deleted 
POLY2  has  not  been  exite 
POLY3  is  entered  at  PT3 

PT4 

POLY3 

POLY  2 

POLY3 

' 

is  exited  and  dele 

•• 

PT5 

POLY2 ,  POLY4 

POLY4 

POLY2 
POLY  4 

is  exited  and  dele 
is  entered  at  PT5 

PT6 

POLY4 

POLY4 

is  exited. 

Fig.  7.16 
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The  next  step  is  to  project  all  the  members  of  NINT  onto 

i 

NI  -  NJ.  The  algorithm  from  here  on  will  be  explained  by 
example.  Referring  to  Fig.  7.16,  let  the  points  in  3-space 
found  by  projecting  PT1  through  PT6  on  NI  -  NJ  be  Nl  through 
N6  respectively.  Each  pair  of  consecutive  nodes  Nl  and  N2, 

N2  and  N3 ,  and  so  on  till  N5  and  N6,  represents  an  edge.  Each 
edge  NINT^  -  NINTI+^  ,  1  I  4=.  5,  is  checked  for  obscurity 

by  the  face  or  faces  whose  projected  polygons  are  recorded  in 
the  Ith  entry  of  FINT.  We  take  the  first  face  of  this  entry 
and  substitute  the  coordinates  of  NINT^.  and  NINTI  +  ^  in  the 
equation  of  the  plane  continuing  the  face  and  get  VALNI  and 
VALNIPl  respectively.  The  coordinates  of  VP  are  also  sub¬ 
stituted  in  the  equation  to  get  VALVP.  If  all  three  values 
have  the  same  sign,  the  edge  is  not  hidden  by  the  face  in 
question.  If  VALNI  and  VALNIPl  have  the  same  sign,  opposite 
to  that  of  VALVP,  the  edge  is  hidden.  If  one  endpoint  is  on 
the  plane  and  the  other  hidden  from  VP,  the  edge  is  hidden, 
and  if  the  other  is  visible  from  VP,  the  edge  is  not  hidden. 

If  both  endpoints  are  on  the  plane,  the  edge  is  again  not 
hidden.  In  all  other  cases  there  is  a  point  of  intersection 
that  must  be  geometrically  found.  Note  that  this  simple 
testing  procedure  results  from  the  fact  that  each  edge 
NINTj  -  NINTI+-j_,  when  projected  onto  the  viewing  plane,  does 
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not  extend  past  the  boundaries  of  the  polygon  recorded  in 
entry  I  of  FINT.  If  such  an  edge  pierces  a  face,  the 
visible  part  of  the  edge  is  that  part  from  the  intersection 
point  in  3-space  to  the  visible  endpoint.  If  the  Ith  entry 
of  FINT  contains  a  second  face,  only  that  part  of  the  edge 
that  is  visible  after  testing  with  the  first  face  need  be 
considered.  When  two  endpoints  are  finally  found  to  be 
visible,  they  are  checked  for  coincidence,  and  if  they  happen 
to  be  coincident  the  edge  is  discarded.  The  very  last  test 
that  has  to  be  made  is  to  find  the  midpoint  of  the  edge  and 
check  it  for  obscurity  using  the  same  methods  and  shortcuts 
of  Phase  II.  If  this  point  is  hidden,  the  edge  is  discarded. 
This  final  test  weeds  out  situations  such  as  those  in  Fig.  7.17. 

When  an  edge  is  definitely  found  to  be  visible,  the  end¬ 
points  are  passed  to  a  routine  that  creates  instructions  to 
draw  the  edge  on  a  graphics  output  device. 
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According  to  the  algorithm,  PTl  and  PT2  will  be  visible. 
It  is  necessary  to  test  MIDPT  in  order  to  discover  that  the 
edge  is  hidden. 

Fig.  7.17 


After  all  edges  of  the  symbols  of  a  picture  have  been 


processed,  Phase  IV  is  started 
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PHASE  IV 

This  phase  fills  in  what  I  call  the  creases  of  a  picture. 
Definition 

A  crease  is  the  common  edge  NI  -  NJ  of  two  intersecting 

faces  Fk  and  FL,  such  that  if  FK  £  F(ST)  ,  FL  £  F(sg^  '  ST  £ 

Su  £  S,  then  neither  NI  nor  NJ  are  original  nodes  of  ST  or  Sy. 

Every  time  a  visible  piercing  point  is  found  during  the 
course  of  the  algorithm  up  till  now,  an  entry  is  made  in  a 
list  which  we  shall  arbitrarily  call  PIERCELIST.  Every  entry 
contains  three  elements: 

1)  the  on-view  face  being  pierced; 

2)  the  on-view  face(s)  to  which  the  piercing  edge  belongs; 

3)  the  coordinates  of  the  piercing  point. 

Since  an  edge  may  belong  to  two  on-view  faces,  two 
entries  may  be  made  in  2)  for  a  case  of  piercing. 

PIERCELIST  is  processed  using  the  following  rules: 

a)  If  for  any  entry  having  FK  as  the  first  element  and  Fj 

as  a  second  element,  there  is  another  entry  having  F^  as  the  second 
element  and  Fj  as  the  first  element,  the  piercing  points  of  the 
two  entries  are  to  be  joined. 

b)  If  any  two  entries  have  identical  first  and  second  elements 
the  piercing  points  are  to  be  joined. 

In  the  following  examples,  NPIERCE  =  ^NPIERCE-^  ^PIERCE^.. 

. o . ,NPIERCEN  \  where  NPIERCE  is  the  set  of  visible  piercing  points. 
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EXAMPLE  I 

V 


Entries  in  PIERCELIST  for  Fig.  7.18 


PIERCED  FACE 
F3 
F3 
FI 
F2 


PIERCING  FACE 
FI 
F2 
F3 
F3 


PIERCING  POINT 
NPIERCE3 
NPIERCE3 
NPIERCE1 

npierce2 


Using  a)  NPIERCE3  has  to  be  joined  to  NPIERCE-^  and  also  to 
NPIERCE2.  Fig.  7.19  shows  the  result. 


Fig.  7.19 


EXAMPLE  II 
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1 


Fig.  7.20 

Entries  in  PIERCELIST  for  Fig.  7.20 


PIERCED  FACE 

FI 

FI 

FI 

F3 


PIERCING  FACE 

F2 

F2 

F3 

FI 


PIERCING  POINT 
NPIERCE1 

npierce2 

npierce2 

npierce3 


From  b)  NPIERCE-|_  should  be  joined  to  NPIERCE2,  and  from  a)  , 
NPIERCE2  should  be  joined  to  NPIERCE3  (see  Fig.  7.21). 
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Fig.  7.21 

The  algorithm  for  generating  creases  falls  down  in  certain 
cases.  Its  great  advantage  is  that  it  consumes  very  little  time 
since  it  only  makes  use  of  by-products  of  the  prior  processing. 
Fig.  7.22,  7.23,  and  7.24  provide  an  example  of  its  inefficiency. 
At  the  time  this  thesis  is  being  written,  a  more  complete 
algorithm  is  being  sought. 
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fl  'x'f  <-&» 


Fig.  7.22 


The  above  configuration  should  result  in  Fig.  7.23,  but 


the  actual  result  is  Fig.  7.24 
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('  £  C\  c.  £ 


Fig.  7.23 


(1*1  &  (Z<.£ 


Fig.  7.24 


CHAPTER  8 


Suggestions  for  Improvement 
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For  the  sake  of  completeness,  brief  mention  is  here 
made  of  certain  areas  of  the  thesis  that  leave  room  for 
improvement.  Remarks  in  passing  throughout  the  thesis 
description  have  probably  already  brought  the  reader's 
attention  to  these  areas. 

1.  The  provision  of  a  good  interface  between  raw  statistics 
and  the  movie-making  system,  with  allowance  given  for  varied 
functional  relations  between  the  statistics  and  the  size  of 
the  symbols  that  represent  them. 

2.  A  capability  for  handling  several  differently  shaped 
symbols  for  maps  carrying  more  than  one  kind  of  demographic 
information. 

3.  A  capability  for  handling  concave  solids  as  symbols. 

4.  Hidden  line  removal  from  the  map  outline  (including 


regional  boundaries) 


APPENDIX 


Was  I  clever  enough?  Was  I  charming 
Did  I  make  at  least  one  good  pun? 

Was  I  disconcerting?  Disarming? 

Was  I  wise?  Was  I  wan?  Was  I  fun? 

-  John  Updike 
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The  Intersection  of  Two  Line  Segments 


Since  the  test  for  determining  whether  two  line  segments 
intersect  is  used  so  frequently,  a  quick  method,  that  is  here 
described,  is  essential.  The  argument  follows  that  outlined 


m 


12 


Consider  two  line  segments.  The  endpoints  of  one  are 
(XI,  Yl)  and  (X2,  Y2) ,  while  the  endpoints  of  the  other  are 
(Ul,  VI)  and  (U2,  V2)  . 

The  first  requirement  for  finding  the  point  of  inter¬ 
section  of  the  line  segments  is  to  determine  whether  such  a 
point  exists.  Following  a  solution  to  this  problem  in  jjL2j 
we  look  at  their  parametric  representations  and  solve  the 


resulting  equations  simultaneously. 

f~X2  XI 
Y2  Yl 


X 

Y 


where  0  ^  T  <1 


U 

V 


,  U2  Ul 
V2  VI 


L1“T  i 

r  pH 

L1-PJ 


where  0  ^  P  ^  1 


There  is  only  an  intersection  if  there  is  a  Tq  and  PQ 
such  that  0  ^  (T0,P0)^  1  ,  resulting  in. 


X 

= 

U 

Y 

V 
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Thus , 


— 

r~ 

—i  -1 

- L 

— I 

To 

s  X2 

-  XI 

U1 

-  U2  1 

U1 

-  XI I 

Lpoj 

’  Y2 

-  Yl 

VI 

-  V2_j 

VI 

-  Yl  i 

_  J 

if  the  inverse  exists.  This  simplifies  to 

Tq  =  (  (VI  -  V2)  .  (U1  -  XI)  -  (U1  -  U2 )  .  (VI  -  Yl)  )  /D, 

PQ  =  (  (X2  -  XI).  (VI  -  Yl)  -  (Y2  -  Yl).  (Ux  -  X±)  )  /D 

where  D  is  the  determinant 

(X2  -  XI).  (VI  -  V2)  -  (U1  -  U2).  (Y2  -  Yl) 

The  problem  of  finding  whether  the  two  segments  intersect 
is  therefore  as  follows: 

a)  Determine  the  values  '  ' 

a  =  X2  -  XI 
b  =  U1  -  U2 
c  =  Y2  -  Yl 
d  =  VI  -  V2 

If-D  =  ad  -  cb  is  equal  to  0,  the  lines  are  parallel.  Otherwise, 
find 

e  =  U1  -  XI 

f  =  VI  -  Yl 

Then  calculate 

T0  =  (de  -  bf)/D 
PQ  =  (af  -  ce)/D 

If  0  ^  (T0,P0) 1  there  is  a  point  of  intersection. 

If  D  is  0  and  the  lines  are  parallel,  they  can  be  checked 
for  collinearity  by  determining  the  value  of  the  expression 
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( Yl  -  VI)  .  (U2  -  Ul)  -  (XI  -  Ul)  .  (V2  -  VI) 

If  this  turns  out  to  have  the  value  0  then  the  lines  are 
collinear.  However,  they  may  still  not  be  superimposed.  If 
(Xg,  Yg)  and  (X 2 /  Y2)  are  the  coordinates  of  the  end  points 
of  one  of  the  two  line  segments,  and  if  (Ug,  Vg)  and  (U2,  V2) 
are  the  coordinates  of  the  endpoints  of  the  other  line  segment, 
then  superimposit ion  occurs  if 


Lpi 

u  p2  u 

(px  p2) 

J  ^  1  where 

xg 

ui 

-  u2 

-  u2 

if  Ui  ji 

U2  for  i  =  1 ,  2 (  or 

=  Yi 

-  V2 

if  Vg  ^ 

V 2  for  i  =  1,  2 

-  V2 

Given  that  the  line  segments  intersect,  the  coordinates 
(INTX,  INTY)  of  the  point  of  intersection  are  given  by 


INTX  = 

XI  + 

To- 

INTY  = 

Yl  + 

V 
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B .  The  Intersection  of  A  Line  With  A  Plane 

Given  a  plane  with  three  points  on  it  having  coordinates 
(XI,  Yl,  Zl)  ,  (X2,  Y2 ,  Z2)  ,  and  (X3,  Y3,  Z3)  ,  and  also 
given  a  line  with  endpoints  (Ul,  VI,  Wl)  ,  and  (U2,  V2 ,  W2)  , 
then  the  point  of  intersection  may  be  found  as  follows. 

Let  the  coordinates  of  the  point  of  intersection  be 
(X,  Y,  Z) .  Since  this  point  is  on  the  line,  we  can  say: 

X  =  Ul  +  C.L 

Y  =  VI  +  C.M 

Z  =  Wl  +  C.N 

where  L,  M  and  N  are  the  easily  calculated  direction  cosines 
of  the  line,  and  C  is  an  unknown  constant.  We  may  write  these 
equations  as  follows: 


(1) 


It  is  necessary  to  find  the  value  of  C,  and  substitution  in 
(1)  will  furnish  the  coordinates  of  the  point  of  intersection. 
The  equation  of  the  plane  is: 


D 


(2) 
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where 


Y1 

Z1 

1 

A  = 

Y2 

Z2 

1 

Y3 

Z3 

1 

Z1 

XI 

1 

Z2 

X2 

1 

Z3 

X3 

1 

XI 

Y1 

1 

X2 

Y2 

1 

X3 

Y3 

1 

!  XI 
| 

Y1 

Z1 

D  =  j  X2 

Y2 

Z2 

!  X3 

Y3 

Z3 

Solving 

(1)  and  (2) 

simultaneously 

\ 

gives 

C  = 

( D  -  /ui\ 

[A\  / 1 

L  'v- 

(  A  \ 

vi ! 

.  B  /  I 

M  . 

!  B 

V  \wi  y 

W///  \ 

N  / 

\c  ! 

X,  Y, 


and  Z  are  now 


obtained  by  substituting 


the  value  of 


C  in  (1) 
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C.  Perspective  Projection 


and 

from 

L25] 

r~  ~ i  r~  — 

Perspective  projection  is  discussed  in  1 2  3j  ,  j 24 
.  The  model  presented  here  is  largely  taken 


The  system  described  in  this  thesis  has  all  points  in 
3-space  provided  in  terms  of  a  rectangular  coordinate  system. 

The  coordinates  (VX,  VY,  VZ)  of  a  viewpoint  VP  may  be  trans¬ 
formed  into  the  spherical  coordinates  (OVER,  UP,  RAD) ,  where 
OVER  and  UP  are  angles  in  degrees,  and  RAD  is  the  radial 
distance/  in  inches,  of  VP  from  the  origin  (see  Fig.  9.1).  A 
picture  plane  must  be  chosen  perpendicular  to  the  line  of 
sight  from  the  viewpoint  to  the  origin,  and  in  such  a  way 
that  all  points  belonging  to  a  visual  scene  are  on  the  other 
side  of  the  plane  from  VP.  Perspective  projection  simply 
involves  finding  the  points  of  intersection  with  the  picture 
plane  of  the  lines  joining  the  viewpoint  to  those  points 
in  3-space  making  up  the  visual  scene. 

Demonstration  of  projection  formulae  begins  by  considering 
the  line  of  sight  along  the  Z-axis,  from  negative  Z  to  the 
origin  in  our  system.  Let  P (X,  Y,  Z)  be  a  point  to  be  projected 
on  the  viewing  plane.  If  F  and  R  are  vectors  whose  magnitudes 
represent  the  focal  length  and  the  distance  of  VP  from  the  origin. 


108 


respectively,  then  the  coordinates  of  the  projected  point  PP(Xp,Yp)are 

XP  =  X.F  Yp  =  Y. F  (see  Fig.  9.2) 

R-Z  R-Z 

We  now  rotate  the  viewpoint  and  projection  plane  by 
UP  degrees  about  the  X-axis,  and  by  OVER  degrees  about  the 
Y-axis.  The  values  of  Xp  and  Yp  now  become: 

Xp  =  _ F  (X  cos  OVER  -  Z  sin  OVER) _ _ 

R  -  Y  +  sin  UP  +  (cos  UP) (Z  cos  OVER  4-  X  sin  OVER) 

Yp  =  F  Py  cos- up  -  (sin  UP)  (Z  cos  OVER  +  X  sin  OVER)l 

R  -  Y  4-  sin  UP  4-  (cos  UP)  [Z  cos  OVER  4-  X  sin  OVER) 

In  viewing  an  object,  all  its  points  have  to  be  projected 
onto  the  picture  plane  using  these  formulae,  and  if  the  object 
is  centred  at  the  origin,  that  is,  if  the  line  of  sight  is  to 
pass  through  the  middle  of  the  object,  the  result  would  be 
very  similar  to  the  way  the  object  would  actually  be  seen  by 
a  person.  If  the  object  is  far  from  the  origin,  or  if  the 
picture  plane  is  too  close  to  the  viewpoint,  a  certain  amount 
of  distortion  will  occur.  For  our  application,  in  creating 
perspective  projection  the  line  of  sight  is  considered  to 
pass  through  the  middle  of  each  frame.  Symbols  are  projected 
onto  the  picture  plane  after  they  have  been  rotated  in  3-space 
to  the  required  orientation,  and  before  the  process  of  hidden 
line  removal.  Since  this  process  makes  use  of  the  depth  of 
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nodes  in  3-space,  the  original  Z-coordinates  are  retained 
after  perspective  transformation  of  the  X  and  Y  coordinates. 
During  hidden  line  removal,  whenever  a  piercing  point  is 
discovered,  it  must  be  subjected  to  perpsective  transformation 
before  it  can  be  used  in  further  processing. 
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D .  Sample  Output 

The  implemented  form  of  the  dynamic  data  display  system 
was  used  to  make  a  1530  frame  (64  second)  movie  of  the  change 
in  the  average  monthly  count  of  the  air  pollutant  SC>2  at  five 
different  locations  in  Toronto  over  a  17  month  period. 

Selected  frames  from  this  movie  are  shown  in  Figs.  9.3, 

9.4,  9.5  and  9.6.  The  frame  numbers  are  shown  above  or  below 
each  frame.  The  reader  will  notice  that  the  viewpoint  is 
slowly  rotated  in  3-space  so  that  he  moves  around  the  map, 
and  also  gently  rises  over  it  from  an  initial  full  side  view 
Perspective  projection  was  not  used  for  this  movie,  a 
deficiency  that  does  little  to  dampen  the  effect  of  the  resu 

There  are  two  steps  to  creating  a  movie,  namely  (a)  generate 
a  condensed  tape  of  graphic  instructions,  and  (b)  plot  the 
tape  created  by  (a) .  The  time-cost  statistics,  excluding 
charges  for  peripherals,  are: 
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CPU  TIME  COST 

CREATE  CONDENSED  TAPE  OF  01//10//26  $440. 2 ( 

GRAPHIC  INSTRUCTIONS 

PLOT  ABOVE  TAPE  00//14//53  55.7(  = 

TOTAL  COST . ' . $496. 0!3 


The  core  required  for  the  dynamic  display  system  is  80. 6K. 

The  resultant  cost  is  well  within  acceptable  bounds 
considering  the  cost  of  using  traditional  animation  methods. 

It  is  true  that  only  5  symbols  are  being  used,  and  also  that 
facilities  for  labels,  a  time  clock,  titles,  visual  scaling, 
etc.  are  not  incorporated,  but  at  $500/min.  their  inclusion 
would  have  to  raise  the  cost  very  sharply  for  the  expense 
involved  in  using  the  system  to  become  prohibitive. 


Frame  101 


Frame  451  Frame  501  Frame  551 


a  me  .,r  F^ame  850  Frame  900 


7rame  1,100  Frame  1,150  Frame  1,200 


Frame  1 , 2  >0  Frame  1,300  Frame  1,330 
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Three  further  examples  from  different  movies  are  provided 

i 

in  Figs.  9.7,  9.8,  and  9.9. 

Fig.  9.7  demonstrates  hidden  line  removal  for  edges  with 
both  endpoints  hidden  and  edges  with  one  endpoint  hidden  and 
the  other  visible. 

Fig.  9.8  demonstrates  the  particular  case  of  an  edge  with 
both  endpoints  visible,  but  having  part  of  that  edge  hidden 
between  the  endpoints. 

Fig.  9.9  shows  piercing  points  and  creases. 
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Fig .  9.7 
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